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Preface 

This  study  is  an  attempt  to  present  analytical  tech¬ 
niques  which  may  be  used  to  evaluate  handling  qualities  of  an 
aircraft  during  the  landing  and  approach  phases  of  flight  re¬ 
gime  in  terms  of  pilot  opinion,  and  to  identify  the  landing 
approach  as  a  closed-loop  tracking  task.  Three  analytical 
techniques:  open-loop  frequency  response  analysis,  closed- 

loop  analog  simulation,  and  open-loop  impulse  response  analy¬ 
sis  are  presented  and  the  results  thus  obtained  are  described. 
It  is  my  hope  that  this  analytical  study  may  be  of  some  help 
to  future  investigators. 

I  wish  to  express  my  indebtedness  to  Ronald  0.  Anderson 
for  his  sponsorship,  expert  assistance,  and  timely  guidance 
in  all  phases  of  this  study.  I  would  like  to  give  special 
thanks  to  my  thesis  advisors  Professor  John  J.  D'Azzo  and 
Captain  James  T.  Silverthorn  for  their  continuous  encourage¬ 
ment,  enthusiasm,  many  hours  of  useful  discussion  and  advice, 
and  above  all  unlimited  patience  throughout  the  course  of 
this  study.  My  sincere  appreciations  are  due  to  Ma jor  J.G.Reid 
for  his  continued  interest  and  providing  publications  in  the 
field  of  pilot  modelling. 
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2 

second,  feet  per  second  ) 

q,q  Perturbation  in  pitch  rate  and  acceleration  (degrees  per 

2 

second,  feet  per  second  ) 

e, e  Perturbation  in  pitch  angle  and  pitch  rate  (degrees, 
degrees  per  second 

a, a  Perturbation  in  angle  of  attack  and  angle  of  attack  rate 
(degrees,  degrees  per  second) 

uQ  Equalibrium  forward  velocity  (feet  per  second) 

WQ  Equilibrium  downward  velocity  (feet  per  second) 

Equilibrium  pitch  angle  (degrees) 

V,p  True  airspeed  (feet  per  second) 

Xu  Dimensional  variation  of  X-force  with  forward  speed  (second  ^ 

Z  Dimensional  variation  of  Z-force  with  forward  speed  (second  ^ 

M  Dimensional  variation  of  pitching  moment  with  forward  speed 

(feet-  second  ) 

X  Dimensional  variation  of  X-force  with  downward  velocity 
w  _1 

(second  ) 

Z  Dimensional  variation  of  Z-force  with  downward  velocity 

w  -1 

(second  ) 
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M  Dimensional  variation  of  pitching  moment  with  downward 
w  _  i 

velocity  (feet  second  ) 

Mq  Dimensional  variation  of  pitching  moment  with  pitch  rate 
(second-^- ) 

2 

g  Acceleration  due  to  gravity  (feet  per  second  ) 

Xx  Dimensional  variation  of  X-force  with  elevator  angle 

ft  2 

e  (feet  per  second^) 

Z,  Dimensional  variation  of  Z-force  with  elevator  angle 

o  2 

e  (feet  per  second  ; 

Mx  Dimensional  variation  of  pitching  moment  with  pitch  rate 
o  0 

e  (second-^) 

«e  Elevator  deflection  (degrees) 
s  Laplace  operator 

h,h  Perturbation  altitude  and  sink  rate  (feet,  feet  per  second) 
R,R  slant  range  and  range  rate  (feet,  feet  per  second) 
u>c  Crossover  frequency  (radians  per  second) 

Kp  Pilot  gain 

t *  Pilot  reaction  time  delay  (seconds) 
t  Pilot  effective  time  delay  (seconds) 

Pilot  Lead  time  constant  (seconds) 

Tj  Pilot  Lag  time  constant  (seconds) 

Neuromuscular  lag  time  constant  (seconds) 

PR  Pilot  Rating 

2 

ISE  Integral  of  Squared  Error  (Volts  -  seconds) 

Phase  Margin  (Degrees) 

Phase  Angle  (Degrees) 

Slope  of  the  phase  curve  (degrees  per  rad/sec) 
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A  4> 

A  i± 
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C/S 
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deviation  In  Pilot  Rating 
Deviation  in  Phase  Angle  (Degrees) 

Deviation  in  slope  of  the  phase  curve  (degrees  per 
rad/sec) 

Control  System 

Aircraft 

Commanded  Pitch  Angle  (Degrees) 

Commanded  Flight  Path  Angle  (Degrees) 

Error  Voltage  in  Pitch  Angle  (Volts) 

Error  Voltage  in  Flight  Path  Angle  (Volts) 

Pilot  gain  in  pitch  loop 

Pilot  gain  in  flight  path  angle  loop 

Time  for  ten  percent  of  final  value  (seconds) 
Duplicating  Time  (seconds) 

Settling  Time  (seconds) 

Pilot  Induced  Oscillations 
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Abstract 

The  approach  and  landing  longitudinal  flying  qualities 
data  generated  by  Calspan  using  variable  stability  NT-33  air¬ 
craft  combined  with  significant  control  system  dynamfcis  was 
analyzed.  Three  different  approaches  were  used  to  determine 
the  relationship  between  pilot  ratings  and  the  pilot  workload, 
system  sensitivity,  and  the  system  performance. 

An  open-loop  frequency  response  analysis  showed  that  for 
a  particular  cross  over  frequency  and  phase  margin  (fixed 
performance),  the  phase  angle  (pilot  workload)  and  slope  of 
the  phase  curve  (system  sensitivity)  at  the  crossover  frequen¬ 
cy  correlated  well  with  the  pilot  ratings. 

An  optimum  pilot  lead  time  for  pitch  tracking,  flight 
path  angle  tracking,  and  both  pitch  and  flight  path  angle 
tracking  tasks  is  determined  by  a  closed-loop  analog  simula¬ 
tion  using  integral  error  squared  as  a  performance  measure. 

A  correlation  is  found  between  pilot  ratings  and  performance. 
An  attempt  is  made  to  explain  the  landing  task  as  a  closed- 
loop  pitch  and  flight  path  angle  tracking  task. 

The  open-loop  impulse  response  of  the  pitch  attitude 
showed  that  the  initial  delay  and  the  transient  behavior  of 
the  pitch  response  have  a  direct  impact  on  the  pilot  ratings, 
pilot  induced  oscillations,  and  comments. 
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ANALYSIS-  OF  THE  EFFECTS  OF  HIGHER 


ORDER  CONTROL  SYSTEMS  ON  AIRCRAFT  APPROACH 
AND  LANDING  LONGITUDINAL  HANDLING  QUALITIES 


I  Introduction 


Background 

The  human  being  is  a  necessary  and  integral  part  of  many 
closed-loop  control  svstems.  In  the  past  two  decades,  a  large 
amount  of  research  effort  has  been  done  to  determine  the  role 
and  response  of  human  beings  in  man-machine  systems.  Mathema¬ 
tical  models  representing  human  behavior  while  perfoi'ming  a 
particular  task  have  been  developed.  The  United  States  Air 
Force  has  made  significant  contribution  in  this  area  of  re¬ 
search  because  of  the  relationship  between  the  pilot  and  the 
aircraft  he  controls. 

The  Air  Force  Flight  Dynamics  Laboratory  at  Wright- 
Patterson  Air  Force  Base  has  been  conducting  studies  and  con¬ 
tracting  for  research  in  developing  a  model  that  would  predict 
pilot  response  in  various  closed-loop  control  tasks.  Models 
for  tasks  such  as  roll  and  pitch  attitude  control  have  been 
successfullv  developed  and  demonstrated  (Refs  4,5).  Of  par¬ 
ticular  interest  at  present  is  the  study  of  aircraft  handling 
aualities  and  pilot  behavior  during  the  glide  slope,  flare, 
and  landing  phases  of  the  flight  regime  in  the  presence  of 
higher  order  control  svstem  dynamics.  The  evaluation  of 


aircraft  handling  qualities  by  the  pilot  during  in-flight 
tests  can  provide  a  means  for  understanding  pilot  behavior, 
and  for  properly  interpreting  the  pilot  ratings  and  comments. 

The  results  of  such  an  evaluation  may  be  used  as  a  guideline  in 
making  qualitative  assessments  of  handling  qualities  of  any  air¬ 
craft  . 

During  the  development  and  in-flight  testing  of  recent 
sophisticated  fighter  aircraft  it  has  become  more  apparent  that 
the  final  landing  task  is  critical  and  must  be  considered 
during  design  and  flight  test  phases.  The  demand  for  increased 
fighter  capability,  coupled  with  present  day  electronic  systems, 
has  made  the  flight  control  systems  more  complex.  The  addition 
of  significant  control  system  dynamics  can  considerably  change 
the  flying  qualities  of  the  aircraft,  and  may  require  extensive 
concentration  and  workload  by  the  pilot  to  achieve  the  desired 
performance . 

The  Flight  Research  Branch  of  Calspan  Advanced  Technology 
Center  conducted  a  study  of  the  effects  of  control  system 
dynamics  on  fighter  approach  and  landing  longitudinal  flying 
qualities  using  an  NT-33A  variable  stability  aircraft.  This 
study  consisted  of  in-flight  simulations  to  generate  an 
approach  and  landing  longitudinal  flying  qualities  data  base 
using  a  highly  augmented  fighter  aircraft  and  contained  docu¬ 
mented  pilot  comments  and  ratings  on  each  flight  (Ref  2). 

The  aircraft  short  period  damping  ratio  and  natural  fre¬ 
quency  were  varied  to  select  five  basic  aircraft  con¬ 

figurations.  Fourteen  control  system  configurations 
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(3  lead  and  10  lag)  ranging  from  first  order  lead  to  fourth 
order  lag  and  the  bare  airframe  dynamics  (control  systems  =1) 
were  selected  and  used  in  various  combinations  to  form 
aircraft-control  system  dynamics.  The  exact  configurations  are 
listed  in  appendix  A  and  the  de'ails  can  be  found  in  Ref  2. 

The  pilot  comments,  pilot  ratings,  and  summary  of  experimental 
results  of  each  configuration  flown  is  also  documented  in  this 
report.  One  of  the  results  of  the  Calspan  study  shows  that 
the  pitch  attitude  tracking  task  performance  is  representative 
of  the  actual  landing  task  performance. 

McDonnell  Aircraft  Company  used  the  data  generated  by 
Calspan  and  applied  the  concept  of  equivalent  systems,  i.e. 
matching  higher-order  aircraft  responses  with  low-order 
equivalents  (Ref  3).  Also,  a  regression  analysis  was  carried 
out  to  develop  a  pilot  rating  expression.  The  study  shows  that 
the  regression  analysis  of  the  equivalent  lower-order  system 
parameters  gave  adequate  prediction  of  pilot  rating  and  the 
equivalent  system  time  delay  is  a  definite  correlating  para¬ 
meter  with  pilot  rating. 

Earlier  studies  conducted  by  McRuer  and  Ashkenas  reveal 
that  short-period  control  of  attitude  and  altitude  using  the 
elevator  can  be  seen  as  either  a  parallel  or  series  loop 
closure  (Ref  6).  In  the  parallel  loop  closure  the  pitch 
and  altitude  channels  are  closed  separately  and  combined  to¬ 
gether  to  generate  the  required  elevator  command.  In  the 
series  loop  closure  the  altitude  is  fed  back  as  an  outer 
loop  to  generate  pitch  command  and  the  pitch  attitude  is  fed  - 
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back  as  an  inner  loop  to  generate  an  elevator  command.  The 

series  closure  is  found  to  be  beneficial  for  low  w  ,  because 

sp 

an  inner  loop  attitude  lead  is  helpful  in  the  outer  altitude 
loop.  The  series  closure  is  detrimental  for  high  “gpj  because 

the  lag  in  the  attitude  loop  degrades  the  outer  altitude  loop. 

Anderson  postulated  that  Pilot  Rating  is  closely  related 
to  pilot  workload  and  closed-loop  system  performance  (Ref  4). 
Specifically,  he  associated  the  pilot  workload  as  the  amount 
of  lead  the  pilot  must  generate  to  maintain  reasonable  response 
and  he  defined  r.m.s.  error  as  the  performance  criterion. 

Problem 

During  the  glide  slope,  flare  and  landing  phases  of  the 
flight  regime,  both  the  pilot  behavior  as  a  feed-back  element 
in  the  closed-loop  pilot-aircraft  control  system  and  the  effect 
of  aircraft-control  system  dynamics  on  the  handling  qualities 
are  as  yet  undefined.  Fundamental  questions,  based  on 
Anderson’s  postulates,  which  need  to  be  answered  are:  1)  What 
is  the  relationship  between  "pilot  workload"  required  to  con¬ 
trol  a  given  set  of  aircraft-control  system  dynamics  during 
landing  task  and  the  corresponding  pilot  rating?  2)  What  are 
the  aircraft  output  variables  the  pilot  is  primarily  interested 
in  controlling,  and  the  relationship  between  the  "performance" 
the  pilot  can  achieve  while  operating  a  given  set  of  aircraft- 
control  system  dynamics  during  landing  task  and  the  correspon¬ 
ding  pilot  ratings?  3)  What  specific  aircraft  dynamic  behav¬ 
ior  does  the  pilot  like  or  dislike  and  reflect  in  his  comments? 


The  Objectives 

The  purpose  of  this  thesis  is  to  carryout  a  qualitative 
and  quantitative  analysis  of  the  data  generated  by  Calspan 
aimed  at  determing  the  effects  of  control  system  dynamics  on 
aircraft  handling  qualities.  This  analysis  provides  an  under¬ 
standing  of  pilot-aircraft  interaction  during  the  landing  task 
in  the  presence  of  significant  control  system  dynamics  and 
answers  the  three  questions  stated  above. 

The  specific  objectives  of  this  thesis  are  the  following: 

1)  To  analyze  the  different  aircraft  and  control  system 
combinations  in  the  frequency  domain  and  to  find  the  correla¬ 
tion  of  the  frequency  response  data  with  pilot  ratings.  It  is 
thought  that  the  phase  angle  at  particular  frequencies  will 
indicate  the  amount  of  lead  the  pilot  must  generate  (workload) 
and  will  therefore  correlate  with  pilot  rating. 

2)  To  simulate  the  pilot-aircraft  dynamics  together  with 
the  control  systems  on  an  analog  computer  as  a  closed-loop 
pitch/flight  path  angle  tracking  task  with  the  purpose  of  find¬ 
ing  the  correlation  between  closed-loop  performance  and  the 
pilot  ratings. 

3)  To  analyze  different  aircraft  and  control  system  com¬ 
binations  in  the  time  domain,  by  applying  an  impulse  input, 
and  interpreting  the  pitch  response  with  the  corresponding 
pilot  ratings  and  the  pilot  comments. 

Approach 

Three  different  approaches  are  used  for  correlating  pilot 
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ratings  and  interpreting  pilot  comments.  In  the  first 
approach  (Chapter  3)  open-loop  frequency  response  plots  of 
pitch  angle  and  flight  path  angle  are  generated  for  different 
aircraft-control  system  combinations  which  were  flown  and  eval¬ 
uated  (Ref  2).  A  pilot  reaction  time  delay  of  0.3  seconds  is 
used  in  all  cases.  The  phase  angle  at  3  radians  per  second 
(rad/sec)  and  the  slope  of  the  phase  curve  at  3  rad/sec  are 
tabulated  to  show  the  relationship  with  pilot  ratings. 

In  the  second  approach  (Chapter  4)  a  closed-loop  analog 
computer  simulation  of  pilot,  aircraft,  and  control  system  is 
studied  in  two  phases.  In  the  first  phase  the  closed-loop 
pitch  response  of  one  aircraft  configuration  to  a  step  pitch 
command  is  studied.  The  pilot  gain  and  lead  time  are  varied 
to  obtain  the  best  performance  for  different  control  system 
configurations.  The  results  are  plotted  to  show  changes  in 
the  pilot  rating  due  to  changes  in  control  system,  as  a 
function  of  performance.  The  flight  path  angle  (y)  loop  is 
closed  next  for  a  step  'y'  command  without  any  pitch  feedback 
and  the  above  procedure  is  repeated.  In  the  second  phase  the 
closed-loop  pitch  and  flight  path  angle  responses  are  studied 
as  in  phase  one  for  pilot-aircraft  dynamics  without  adding 
control  systems  but  changing  the  aircraft  dynamics.  The  pilot 
gain  and  lead  time  are  varied  to  obtain  the  best  performance. 
The  results  are  plotted  to  show  the  changes  in  pilot  rating  as 
a  function  of  performance.  In  addition,  using  aircraft  number 
one  and  different  lag  control  systems,  the  effect  of:  closing 
only  the  outer  y  loop  for  a  step  y  command  is  studied. 
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In  the  third  approach  (Chapter  5)  the  open-loop  time  re¬ 
sponses  of  pitch  attitude  to  an  impulse  are  generated.  Some 
of  the  pilot  comments  and  the  pilot  ratings  are  correlated 
with  the  pitch  behavior  in  the  initial  few  tenths  of  a  second. 

For  the  above  analyses,  aircraft  equations  of  motion  and 
a  pilot  describing  function  are  required.  Chapter  II  develops 
these  equations  and  presents  the  concepts  of  pilot  ratings 
applied  towards  this  study.  Anderson's  concept  is  modified  to 
include  "system  sensitivity"  as  a  parameter  in  evaluating  pilot 
ratings,  measured  in  terms  of  the  slope  of  the  phase  angle 
about  the  cross-over  frequency.  A  description  of  the  landing 
task  as  seen  by  a  pilot  is  also  presented. 

Scope 

It  is  emphasized  that  this  thesis  is  neither  an  attempt  to 
develop  a  new  pilot  model  nor  to  develop  an  expression  for 
pilot  rating  which  can  be  readily  used  for  any  given  situation. 
The  study  is  limited  to  the  anlysis  of  the  Calspan  data  to  dem¬ 
onstrate  the  effect  on  pilot  ratings  for  the  pilot-aircraft 
dynamics  as  control  systems  are  varied.  A  "Lead  only"  form  of 
McRuer  model  for  the  pilot  is  used  with  a  0.3  second  reaction 
time  delay.  For  the  analog  simulation  a  short  period  airframe 
model  is  used  and  the  high  frequency  feel  and  actuator  dynamics 
are  neglected.  The  time  integral  of  error  squared  is  used  as 
the  performance  measure.  For  the  open  loop  pitch  attitude  time 
responses  the  aircraft  short  period  dynamics  are  used.  No  wind 
gust  or  random  inputs  are  considered  and  the  aircraft  forward 
speed  throughout  the  landing  task  is  assumed  to  be  constant. 
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II.  Theoretical  Development 


This  chapter  presents  the  theoretical  concepts  and  the 
associated  mathematical  equations  which  are  used  for  the 
analysis  in  the  subsequent  chapters  of  this  report.  The 
linearized  aircraft  longitudinal  equations  of  motion  are 
developed  and  expressed  in  the  Laplace  domain.  The  short 
period  approximation  is  used  to  develop  the  short  period 
aircraft  equations.  A  description  of  the  landing  task,  as 
seen  by  a  pilot,  is  presented  and  the  significance  of  dyna¬ 
mic  behavior  of  the  pitch  attitude,  0,  and  the  flight  path 
angle,  y  ,  during  the  landing  task  is  discussed.  The  pilot 
describing  function  for  a  single  input-single  output  closed- 
loop  tracking  task  is  described  and  the  model  used  for  the 
analog  simulation  is  developed.  The  pilot  rating  concepts 
are  enumerated  and  the  factors  such  as  performance,  work¬ 
load,  and  sensitivity  which  may  be  taken  into  account  by  a 
pilot  while  providing  a  rating  are  explained.  The  methods 
used  for  analysing  the  aii craft-control  system  dynamics  for 
measuring  the  performance,  workload, and  sensitivity  are 
briefly  mentioned. 

The  linearized  aircraft  longitudinal  equations  of  mo¬ 
tion  expressed  in  body  axes,  in  the  absence  of  external 
disturbances,  are:  (Ref  2:211) 
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where : 

u  =  perturbed  forward  speed 
w  =  perturbed  downward  velocity 
q  =  perturbed  pitch  rate 
e  =  perturbed  pitch  angle 
uq  =  Equilibrium  forward  velocity 

Wq  =  Equilibrium  downward  velocity 

eQ  =  Equilibrium  pitch  angle 

For  small  angles,  uq  =  V,p  >  the  true  trim  airspeed,  and 
the  perturbation  angle  of  attack  a  =  which  implies 

• 

*  w  * 

a  =  —  .  Substituting  for  w  and  w  ,  the  perturbation  equa- 

o 

tions  become: 
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Taking  the  Laplace  Transform  of  equations  (2)  with  zero  in¬ 
itial  conditions,  and  expressing  the  resulting  equations  in 
Matrix  format  yields: 
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From  equation  (3)  the  longitudinal  transfer  functions  repre¬ 
senting  perturbations  in  u,a,q  and  e  can  be  derived  for  the 
elevator  input,  6e,as  shown  in  appendix  B. 

The  denominator  of  the  longitudinal  transfer  function 
derived  from  equation  (3)  is  a  fourth  order  polynomial.  In 
general,  the  polynomial  factors  into  two  complex  pair  of 
roots,  one  with  a  relatively  high  damping  ratio  and  natural 
frequency  called  the  short  period  poles,  and  the  other  with 
a  very  low  damping  ratio  and  natural  frequency  called  the 
phugoid  poles.  The  short  period  mode  is  primarily  made  up 
of  perturbations  in  a  and  o  with  very  small  variations  in 


forward  speed.  The  phugoid  mode  consists  of  variations  in 
0  and  u  with  almost  a  constant  <*.  The  phugoid  mode,  due  to 
its  low  damping  and  low  frequency  of  oscillation,  appears 
for  a  longer  period  of  time.  But  pilots  generally  control 
these  oscillations  without  any  difficulty.  During  the  land¬ 
ing  task,  the  perturbation  in  the  forward  speed  is  controll¬ 
able,  either  manually  or  by  an  automatic  speed  control  system. 
Also,  the  time  frame  in  which  the  landing  task  is  completed, 
is  relatively  small.  Therefore,  neglecting  the  phugoid  mode 
during  the  landing  task  will  not  cause  any  serious  deviations 
in  the  overall  results  of  the  analysis  of  the  task.  It  is, 
therefore,  assumed  that  the  small  perturbations  in  the  forward 
speed  are  negligible. 

If  the  assumption  of  constant  speed  and  0Q  =  o  ,  a  rea¬ 
sonable  assumption  for  the  landing  task,  are  made  the  u-degree 

c 

of  freedom  becomes  superfluous  and,  with  the  substitution  of 
q  =  e  ,  equation  (3)  reduces  to:  (Ref  7:6-30) 


Equation  (4)  is  refered  to  as  the  short  period  approximation. 
The  longitudinal  transfer  functions  representing  perturbations 
in  a  and  0  for  an  elevator  input  can  be  derived  from  equation 
(4)  as  shown  in  appendix  B.  The  values  of  the  stability 
derivatives  for  the  five  aircraft  configurations  are  listed  in 
appendix  B. 


The  Landing  Task  As  Seen  By  A  Pilot 

During  the  approach-to-landing  phase  of  flight  the  air¬ 
craft  is  initially  on  the  glide-slope  and  the  pilot  must  de¬ 
cide  on  the  minimum  height  to  flare.  At  this  time  the  pitch 
attitude, e,  and  the  flight  path  angle, y,  appear  to  be  of  con¬ 
siderable  importance  to  the  pilot.  The  minimum  height  to 
flare  is  an  imaginary  window  on  the  glide  slope  which  the 
pilot  selects  relative  to  the  runway  threshold  or  touchdown 
point.  The  pilot  starts  the  flare  maneuver  in  this  window 
in  order  to  achieve  a  smooth  exponential  landing  trajectory. 
The  selection  of  the  imaginary  window  in  space  depends  on 
the  aircraft  forward  speed,  aircraft  dynamic  characteristics, 
and  the  experience  of  the  pilot  in  handling  the  aircraft. 
During  the  flare  maneuver  the  sink  rate,  h  (vertical  velocity) 
and  the  slant  range  rate,  R  (the  rate  at  which  the  runway 
threshold  or  touchdown  point  appears  to  be  approaching)  are 

the  motion  cues  that  the  pilot  gets  through  his  peripheral 

• 

vision.  The  pilot  pitches  up  and  establishes  a  desired  h 

•  • 

and/or  R  until  touch  down,  where  h,h  and  R  go  to  zero.  All 
three  output  variables,  h,  h,  and  R  are  mathematically  related 
to  the  flight  path  angle  y  as  shown  below: 

For  small  angles,  assuming  constant  forward  speed  through 
out  the  landing  task,  the  flight  path  angle  y  is  given  by: 

(Ref  8:15,83) 


Y 


6  —  a 


(5) 


The  sink  rate  h  is  given  by: 

h  =  uQSin  y  =  uqy  (6) 

•  • 

The  slant  range  rate  R  is  the  vector  sum  of  uQ  and  h.  Since 

•  • 

uQis  assumed  constant,  R  is  directly  proportional  to  h  or  t- 

Therefore,  all  the  motion  cues  are  direct  functions  of  flight 
path  angle  y.  Hence  the  dynamic  behavior  of  0  and  y  for  eleva¬ 
tor  commands  (  se)  is  representative  of  the  over  all  landing 

task.  The  aircraft  transfer  functions  e/<5  G  ,  a/s  e  and  y /s  e 

for  all  five  aircraft  configurations  are  developed  in  appen¬ 
dix  B. 

Pilot  Describing  Function 

The  determination  of  a  suitable  pilot  model  is  not  a 
straight  forward  process.  Understanding  human  behavior  and 
his  reactions  in  a  given  situation  has  been  a  serious  con¬ 
cern  in  the  studies  carried  out  so  far.  A  human  being  per¬ 
ceives  his  task,  evaluates  it  mentally,  then  takes  certain 
action  which  is  delayed  to  varying  degree  depending  upon  his 
psychophvsiological  state  and  firmness  of  decision.  This  is 
generally  refered  to  as  a  reaction  time  delay  which  accounts 
for  the  delay  from  perception  to  the  decision  for  corrective 
action.  The  execution  of  his  decision  is  performed  through 
the  neuromuscular  system,  which  not  only  contributes  some 
delay  but  also  high  frequency  attenuation,  usually  approxi¬ 
mated  by  a  first  order  lag.  In  a  single  loop  compensatory 
tracking  task  the  human  being  interfaces  with  the  system  and 


adjusts  himself  to  achieve  the  desired  performance,  just  like 
a  compensator  in  a  control  system.  The  human  compensation 
is  refered  to  as  equalization.  The  equalization  is  repre¬ 
sented  in  terms  of  a  lead/lag  transfer  function.  The  human 
being  is  also  able  to  adjust  the  overall  system  gain  to 
select  a  desired  operating  frequency  and  phase  margin.  The 
operating  frequency  is  known  as  the  cross-over  frequency,  <<>c 

at  which  the  overall  gain  is  adjusted  for  Odb  on  the  open 
loop  frequency  response  plots  for  a  desired  phase  margin. 

The  mathematical  expression  for  the  aforesaid  human 
behavior  in  a  single  loop  compensatory  tracking  task  is 
given  by: 

a 

K  e“T  s  (Tt  S  +  1) 

Y  (S)  =  - 2 - L -  (7) 

P  (TjS  +  1)(TnS  +  1) 

which  is  known  as  the  McRuer  model  (Ref  1:29;  10:8). 

In  this  model  the  pure  time  delay  is  represented  by 

e-T  S  ,  the  gain  by  and  the  general  lead  and  lag  terms  as 

T^  and  Tj  respectively.  The  first  order  lag  time  constant 

approximation  of  the  neuromuscular  system  is  represented  by 
T^.  Often  the  neuromuscular  lag  term  is  treated  as  an  addi- 

^  M  T  S 

tional  delay,  and  the  pure  delay  term,  t  is  modified  as  e  , 

A 

where  T  =  T  +  tn  '  Hence,  the  Laplace  transform  of  the 
simplified  version  of  the  pilot  describing  function  is: 


(8) 


Yp(S> 


V"'S(tls  + 

(TjS  +  1) 


The  model  parameters  t,T^,Tj  and  are  adjusted  to  pro¬ 
duce  both  stability  and  good  closed-loop  performance  without 
excessive  workload.  In  particular  by  increasing  and  T^, 

while  reducing  t,  the  pilot  can  achieve  the  desired  level  of 
stability  and  performance.  This  is  done  at  the  expense  of 
pilot  activity  and  workload.  Furthermore,  the  pilot  can  in¬ 
troduce  lag  (Tj)  to  maintain  desirable  low  frequency  response. 

Obviously  the  adjustment  of  these  parameters  is  very  dependent 
upon  the  system  the  pilot  is  trying  to  control.  One  way  of 
optimizing  performance  is  to  adjust  the  model  parameters  to 
minimize  the  mean  squared  error  as  the  performance  measure¬ 
ment.  A  performance  measurement  is  an  evaluation  of  how  well 
a  system  operates,  how  closely  the  output  follows  the  input, 
or  how  small  the  error  can  be  maintained  over  a  given  time 
period.  The  time  integral  of  error  squared  has  been  success¬ 
fully  used  to  evaluate  the  performance  of  pilot-aircraft  con¬ 
trol  systems  (Ref: 9).  The  details  of  adjustment  rules  for 
the  human  operator's  describing  function  can  be  found  in  the 
literature  on  the  subject  (Ref:  1:30;  10:9). 

For  the  purpose  of  this  thesis  the  pilot  model  used  for 
analog  simulation  is  of  the  form:  (Ref  4) 

Yp(S)  =  Kpe"TS(TLS  +  X) 


(9) 


This  ’’lead  only"  form  is  selected  because  of  the  fact  that 
pilot  lag  may  be  neglected  when  higher  order  effects,  such 
as  control  system  lags,  are  considered.  Besides,  the  short 
period  approximation  of  aircraft  dynamics  provides  a  reason¬ 
ably  good  low  frequency  response  due  to  the  "free  s"  in  the 
denominator  of  the  °(s)/fi  (s)  transfer  function.  Also  it  is 

only  when  the  natural  frequency  of  the  controlled  element  is 
very  large  >  6  rad/sec)  that  the  pilot  lag  may  be  desir¬ 

able  (Ref  4:4).  Since  the  aircraft  under  study  have  short 
period  frequencies  varying  between  1  rad/sec  and  4  rad/sec, 
the  form  of  pilot  describing  function  as  shown  in  equation 
(9)  is  considered  to  be  a  satisfactory  model  for  analog  sim¬ 
ulation.  A  value  of  t  =  0.3  sec  is  used  for  this  study 
(Ref  2). 

Pilot  Rating  Concepts 

Pilot  Ratings  are  based  on  numerical  rating  scales,  such 
as  the  Cooper-Harper  scale,  which  represent  an  attempt  to  re¬ 
late  pilot  comments  about  the  ease  or  difficulty  with  which 
airplanes  can  be  controlled  in  certain  flight  situations.  Un¬ 
fortunately  the  ratings  are  ordinal  scales  subjectively 
applied  and  hence  difficult  to  quantify.  Predictions  of 
pilot  ratings  could  be  made  if  mathematical  relations  were 
developed  between  the  numerical  Tating  scale  and  how  hard  a 
"pilot  must  work"  to  achieve  "desired  closed-loop  performance" - 
Developing  mathematical  relations  appears  to  be  extremely  com¬ 
plex,  though  certain  trends  have  been  identified.  For  example 


a  pilot  objects  seriously  if  he  has  to  generate  "leads"  (T^) 

of  more  than  0.5  to  1.0  seconds  (Ref  7:12-9). 

In  order  to  develop  closed-loop  analysis  procedures  for 
flying  qualities  assessments,  some  generalized  criteria  are 
required  in  terms  of  pilot-aircraft  dynamics,  task  perfor¬ 
mance,  and  the  pilot  workload.  Therefore  the  factors  which 
may  be  taken  into  account  by  a  skilled  test  pilot  in  provi¬ 
ding  an  opinion  or  rating  are:  (Ref  1:35) 

a)  Measure  of  task  performance 

b)  Pilot  workload 

c)  System  sensitivity 

These  factors  obviously  are  affected  by  the  particular 
aircraft  dynamics,  including  any  additional  control  systems. 
As  a  mathematical  function,  pilot  ratings  (PR)  can  be  ex¬ 
pressed  as: 

PR  =  f (Performance ,  workload,  sensitivity,...  )  (10) 

Performance 

Task  performance  is  a  measure  of  how  well  the  system 
responds  to  input  commands.  The  closed-loop  system  quanti¬ 
ties,  such  as  the  response  time  and  the  integral  error 
squared  (ISE),  are  indicative  of  the  accuracy  of  the  task. 

In  general,  both  the  response  time  and  ISE  are  functions  of 
the  cross-over  frequency  (w  )  and  the  phase  margin  (<t>  ). 

c  m 

Ratings  are  known  to  be  poor  if  u>c  and/or  ISE  are  inadequate. 
A  pilot  is  an  element  in  the  closed-loop  system  and  tries  to 
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adjust  w  and  to  maintain  sufficient  phase  margin  to  minimize 
ISE.  Pilots  generally  consider  1  rad/sec  <  <*>c  <  5  rad/sec  and 

40°  <  4>m  <  60°  to  be  satisfactory  operating  regions  for  ad¬ 
equate  performance  (Ref  6:81;  7 • 12 .4-12 .6)  .  The  choice  of 

w  and  <t>  ,  however  depends  on  the  aircraft  control  system 
cm  ' 

dynamics . 

Pilot  Workload 

The  pilot  workoad  can  be  divided  into  pilot  activity 
and  pilot  equalization.  Pilot  activity  is  dependent  on  the 
pilot  gain  Kp  and  is  an  important  component  of  Pilot  Rating 

in  the  presence  of  gust  disturbances  and  remnant.  Since 
gust  disturbances  and  remnant  effects  are  not  considered  in 
this  study,  the  pilot  activity  is  neglected.  Besides  the 
pilot  gain  is  also  determined  by  the  mechanical  arrangement 
of  the  controls  placed  between  the  pilot  and  the  elevator, 
like  the  gearing  ratio.  The  data  generated  by  Calspan  had 
the  facility  for  the  pilot  to  select  the  gearing  ratio  for 
each  flight  evaluation  of  a  configuration  (Ref:2). 

The  pilot  equalization  can  be  represented  in  two  ways 
(Ref  1).  One  way  is  to  determine  the  slope  of  the  pilot's 
amplitude  ratio  evaluated  near  the  cross-over  frequency. 

The  other  way  is  to  use  the  pilot's  "lead"  time  constant.  In 
this  study  the  pilot  "lead"  time  constant  is  used  to  represent 
workload,  as  the  pilot  gain  is  an  adjustable  parameter. 
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Sensitivity 

System  sensitivity  is  a  measure  of  performance  degrad¬ 
ation  due  to  small  deviations  in  the  closed-loop  operating 
frequency  or  <»>  .  It  may  neither  be  possible  nor  desirable 

for  a  pilot  to  maintain  a  fixed  all  the  time  during  the 

execution  of  a  particular  task.  Small  deviations  about  u>c 

cause  variations  in  the  overall  system  phase  angle  or  phase 
margin  as  seen  in  the  frequency  domain.  A  sharp  variation 
of  phase  angle  about  «c  (measured  in  degrees  per  rad/sec) 

will  result  in  considerable  change  in  <t>^,  and  hence  a  signi¬ 
ficant  closed-loop  performance  degradation  in  terms  of  ISE, 
whereas  a  gradual  variation  of  phase  angle  about  ojc  will  not 
degrade  the  system  performance  rapidly.  Therefore,  the  slope 
of  the  phase  angle  curve  about  wc  indicates  the  relative  sy¬ 
stem  sensitivity  to  small  changes  in  t»c. 

Methods  For  Longitudinal  Handling  Qualities  Analysis 

The  mathematical  equations  representing  the  aircraft 
characteristics  are  known,  as  the  formulation  of  transfer 
functions  is  usually  accomplished  during  the  aircraft  design 
phase.  The  transfer  functions  of  additional  control  systems 
designed  to  enhance  the  aircraft  capabilities  can  also  be 
formulated.  Therefore  an  open-loop  frequency  response  plot 
of  the  complete  system  can  be  generated.  The  phase  angle  and 
the  slope  of  the  phase  curve  about  a  nominal  cross-over  fre¬ 
quency  where  pilots  are  known  to  operate  is  a  measure  of  the 
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pilot  workload  and  system  sensitivity.  Closed-loop  pilot- 
aircraft  dynamics  can  be  simulated  on  an  analog  computer  for 
a  given  tracking  task,  assuming  a  form  for  the  pilot  model. 

The  pilot  parameters  can  be  adjusted  to  minimize  the  integral 
of  error  squared,  which  is  used  as  a  performance  measure. 

The  closed-loop  operating  frequency  and  the  phase  margin  thus 
obtained  is  an  indication  of  system  performance,  and  the  pilot 
parameters  are  an  indication  of  estimated  workload.  An  open- 
loop  pitch  attitude  response  to  an  impulse  can  be  generated 
using  the  aircraft-control  system  dynamics.  The  system  be¬ 
havior  in  the  initial  few-tenths  of  a  second  reveals  the  over¬ 
all  system  lag  or  delay  and  can  be  used  to  interpret  the  ex¬ 
pected  pilot  ratings  and  comments.  These  three  approaches: 

open-loop  frequency  response,  closed-loop  analog  simulation, 

* 

and  open-loop  impulse  response  are  analyzed  in  the  next  three 


chapters . 


Ill  Open  Loop  Frequency  Response  Analysis 


The  open  loop  frequency  response  analysis  provides  an 
overall  picture  of  magnitude  and  phase  relationship  of  the  out¬ 
put  of  a  system  with  respect  to  a  sinusoidal  input.  From  the 
frequency  response  plots  the  gain  margin  and  the  phase  margin 
of  the  system  can  be  measured  at  any  desired  operating  fre¬ 
quency.  The  gain  and  the  phase  margins  are  related  to  the 
effective  damping  ratio  of  the  resulting  closed  loop  system, 
and  hence  the  overall  system  transient  behavior  (Ref  11:288). 
The  gain  and  phase  margins  also  provide  a  basis  for  selection 
of  a  suitable  cascade  compensator  in  order  to  achieve  a  desired 
closed-loop  system  performance  (Ref  11:386). 

This  chapter  presents  a  brief  description  of  the  overall 
system.  The  digital  computer  program  "TOTAL"  is  used  for  the 
frequency  response  analysis  (Ref  12)  and  the  results  obtained 
from  the  analysis  are  explained.  Pilot  ratings  are  correlated 
with  the  phase  angle  data  using  least  square  curve  fit.  The 
justification  for  this  analysis  is  derived  from  the  fact  that 
the  desired  phase  margin  at  the  crossover  frequency  for  a  rea¬ 
sonable  closed-loop  system  performance  is  achieved  by  the 
pilot  lead.  Hence  the  pilot  workload  is  directly  related  to 
the  resultant  system  phase  angle  at  the  crossover  frequency 
in  the  frequency  domain.  The  slope  of  the  phase  curve  at  the 
crossover  frequency  is  a  measure  of  the  system  sensitivity. 
Therefore  a  correlation  can  exist  between  pilot  ratings  and 
the  phase  angle  data. 
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Two  methods  of  correlating  pilot  rating  expressions  are  presen¬ 
ted.  In  the  first  method  an  expression  is  developed  to  fit  the 
total  pilot  rating  with  phase  angle.  In  the  second  method  the 
deviation  in  the  pilot  rating  with  respect  to  the  best  rating 
is  correlated  with  the  corresponding  deviations  in  the  phase 
angle  data. 

The  block  diagram  of  the  open-loop  system  under  investi¬ 
gation  is  shown  in  Fig.  1.  A  pilot  reaction  time  delay  of 
0.3  seconds  is  included.  In  the  Laplace  transform  domain,  a 

__  S  T 

pure  delay  of  t  seconds  becomes  a  multiplying  factor  of  e 

Since  |e‘’^“T|  is  unity,  this  does  not  affect  the  gain  plot, 
but  on  the  phase  curve  it  introduces  an  additional  phase  lag 
of  wt  radians  at  each  frequency  w. 


Fig  1  Overall  Open-Loop  System  Block  Diagram 

The  variables  in  the  analysis  are  the  elements  of  the  con¬ 
trol  system  dynamics.  The  aircraft-control  system  configura¬ 
tion  used  and  the  corresponding  pilot  ratings  are  shown  in 
appendix  A.  The  feel  system  and  actuator  dynamics  are  held 
constant  for  all  configurations.  The  transfer  functions  used 
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for  the  control  systems,  feel  system  and  actuator  dynamics 

are  also  shown  in  appendix  A.  Using  the  g—  transfer  function 

e 

and  transfer  function  developed  in  Appendix  B  for  each  air- 

e 

craft  configuration,  the  control  systems  were  varied  to  gener¬ 
ate  open-loop  frequency  response  curves  using  the  digital 
computer  program  "TOTAL".  The  magnitude  and  phase  curves  for 
selected  configurations  are  attached  as  Appendix  E. 

A  close  examination  of  the  phase  angle  plots  shows  some 
distinct  correlating  trends  between  the  pilot  ratings,  and  the 
phase  angle  and  slope  of  the  phase  curve  in  the  vicinity  of 
1  to  4  rads /sec.  Control  systems  generating  more  phase  lag  in 
this  frequency  region  have  a  much  higher  rating  as  compared  to 
the  control  systems  generating  less  lag  in  this  region.  There¬ 
fore,  the  phase  angles  of  each  aircraft-control  system  config¬ 
uration  were  measured  at  1,2,3,  and  4  rads/sec  and  tabulated. 
Tables  El  thru  E  10  (Appendix  E)  show  the  values  of  phase 
angles  $  and  the  change  in  4>  per  rad/sec. 

In  chapter  II  the  pilot  rating  is  explained  as  a  func¬ 
tion  of  performance,  workload,  and  sensitivity.  Since  no  in¬ 
formation  about  any  of  the  three  variables  is  known,  a  fre¬ 
quency  of  3  rad/sec (representing  a  fixed  level  of  performance) 
was  selected  as  the  cross  over  frequency  for  the  purpose  of 
analysis.  The  choice  of  3  rad/sec  is  based  on  the  earlier 
finding  that  the  pilots  like  performing  closed-loop  tasks 
around  3  to  4  rad/sec.  In  order  to  control  the  system  for 
satisfactory  operation,  the  pilot  must  generate  sufficient 
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lead  to  achieve  the  required  phase  margin  at  3  rad/sec.  The 
phase  margin  and  the  cross  over  frequency  are  directly  related 
to  the  performance.  Therefore,  this  analysis  is  based  on  the 
assumption  that,  for  a  fixed  performance,  i.e.,  for  given  o>c 

and  the  phase  angle  and  slope  of  the  phase  curve  at  uc 

are  a  relative  measure  of  pilot  workload  and  system  sensitiv¬ 
ity,  respectively. 

Determination  of  Pilot  Rating  Expression 

It  is  possible  to  correlate  the  pilot  ratings  with  the 
phase  angle  data  recorded  in  tables  El  thru  E10.  Since  the 
aircraft  dynamics  are  varied  to  generate  five  different  con¬ 
figurations,  the  data  is  analysed  to  determine  a  rating  ex¬ 
pression  for  each  configuration  as  the  control  systems  are 
varied.  Anderson's  concept  (Ref  4)  that  pilot  rating  is 
based  on  some  linear  combination  of  performance  and  workload 
is  modified  to  include  the  system  sensitivity  as  a  third 
variable.  Therefore,  for  fixed  performance,  the  pilot  rating 
is  a  linear  combination  of  workload  and  sensitivity.  This 
leads  to  a  rating  expression  of  the  form 

PR  =  a.,<t>-5  +  ,  /ii\ 

1  3  2  6  w  Avg  (11) 

where  and  02  are  some  constants  to  be  determined,  ^  is 

the  phase  angle  of  the  open-loop  system  shown  in  Fig.  1  at 
3  rad/sec  and  yy  |  is  the  average  slope  of  the  phase 

curve  between  2  to  4  rad/sec.  This  average  is  obtained  by 
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taking  the  mean  phase  difference  between  2  and  3  rad/sec, 
and  3  and  4  rad/sec.  The  estimated  pilot  rating  is  represen¬ 
ted  by  PR.  The  procedure  of  linear  least  square  curve  fit  is 
used  to  evaluate  and  such  that  the  square  of  the  error 

between  the  estimated  and  actual  Pilot  Ratings  is  minimized. 
This  procedure  is  discussed  and  the  equations  used  are  de¬ 
rived  in  Appendix  C. 

The  values  of  and  «2  are  obtained  using  the  equations 
developed  in  Appendix  C  and  the  pitch  angle  frequency  response 
data  tabulated  in  Tables  El  thru  E5.  The  data  for  aircraft 
number  3,  listed  in  Table  E3  is  not  used  as  it  has  a  very 
significant  slope  variation  between  3  and  4  rad/sec .  There¬ 
fore,  averaging  the  two  slopes  is  not  a  reasonable  approxima¬ 
tion  for  this  set  of  data.  Also  the  data  obtained  by  adding 
lead  control  systems  to  the  five  aircraft  were  not  included 
because  the  lead  control  system  improves  the  phase  angle  and 
slope,  hence  the  expression  shown  in  equation  11  would  yield 
better  pilot  ratings.  But  the  actual  pilot  ratings  given  by 
the  pilot  for  all  the  aircraft  flown  with  lead  control  systems 
are  always  worse  than  the  aircraft  with  no  control  system  add¬ 
ed.  This  reverse  trend  in  pilot  ratings  could  not  be  explained 
at  this  point  and  therefore  these  data  points  were  omitted. 

The  values  of  and  02  for  aircraft  1,2,4  and  5  obtained 

by  the  linear  least  square  curve  fit  analysis  are  listed  in 
Table  Ell,  Appendix  E.  The  estimated  PR  and  the  actual  ratings 
are  also  tabulated.  Figure  2  shows  a  plot  of  estimated  PR 
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ESTIMATED  PR 


versus  the  actual  PR.  Fourteen  out  of  the  thirty  points 
plotted  are  outside  the  limits  of  ±1  from  the  best  line  of  co¬ 
relation.  Out  of  the  fourteen  points  six  data  points  repre¬ 
sent  aircraft  2,  five  data  points  represent  aircraft  4,  and 
three  represent  aircraft  1. 

In  the  above  analysis,  since  the  average  slope  did  not 
accurately  represent  the  actual  slope  of  the  phase  curve  at 
three  rad/sec,  an  alternate  expression  is  developed  to  in¬ 
clude  both  the  slopes  from  2  to  3  rad/sec  •—  J  2.3  an<*  3  to  4 

rad/sec  4^  ,  /  as  shown: 

610  j-4 

PR  =  a14>3  +  a2  6^  |  2-3  +  a3  6^  |  3-4  (12) 

where  a1>0l2  anc*  a3  are  the  constants  evaluated  applying  the 

method  of  least  squares.  The  data  for  aircraft  3  is  also  in¬ 
cluded.  The  results  thus  obtained  are  tabulated  in  Table  E12. 

A  plot  of  the  estimated  against  the  actual  PR  is  shown  in 
Fig.  3.  From  the  figure  it  is  seen  that  eleven  out  of  thirty 
five  points  are  outside  the  ±1  PR  boundary.  Seven  out  of  the 
eleven  points  represent  aircraft  2,  three  data  points  represent 
aircraft  1  and  one  data  point  represents  aircraft  4.  The 
second  PR  expression  (eqn. 12)  did  show  a  better  corelation  as 
compared  to  the  first  (eqn.  11),  but  of  course  it  added  an  ad¬ 
ditional  degree  of  freedom. 

As  an  alternate  approach  for  correlating  the  phase  angle 
data,  the  deviations  in  the  pilot  ratings  with  respect  to  the 
best  rating  in  the  data  was  curve  fit  to  corresponding 
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Figure  3.  Estimated  Verses  Actual  PR  Using  Pitch  Phase  Data 
And  Equation  12. 


deviations  in  the  phase  angles  and  the  slopes.  Once  again  a 
frequency  of  3  rad/sec  was  selected.  The  PR  expression  thus 
used  was: 


APR 


(  A  <t>  ) 


“2 


Avg 


(13) 


where  and  are  the  unknown  constants,  PR  is  the  devia¬ 
tion  in  pilot  rating  with  respect  to  the  best  PR,  A <t>  is  the 


corresponding  deviation  in  the  phase  angles,  and  A 


6  to 


Avg 


is  the  average  deviation  in  the  slope  with  respect  to  the 
average  slope  corresponding  to  best  rating. 

The  values  of  and  for  each  aircraft,  obtained  by 


using  the  method  of  least  squares,  and  the  estimated  PR  are 
tabulated  in  Table  E13.  A  plot  of  actual  versus  the  estimated 
APR  is  shown  in  Fig.  4.  Eight  out  of  thirty  points  are  outside 
the  ±1PR  region.  Out  of  the  eight  data  points  four  represent 
aircraft  2,  two  represent  aircraft  1,  and  two  represent  aircraft 
4.  This  analysis  primarily  indicates  the  effects  of  lag  control 
systems  on  the  pilot  ratings. 

The  APR  expression  of  eqn.  13,  was  slightly  modified  to 
include  a  bias  (a0)  as  shown  and  the  procedure  was  repeated. 


APR 


a  + 
o 


a^Aifr) 


«2 


Avg 


(14) 


Only  aircraft  1  and  2  were  analyzed  to  see  if  any  noticeable 
improvement  in  the  correlation  exists.  The  results  thus 
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ESTIMATED  APR 


And  Equation  13 


obtained  are  tabulated  in  Table  E14.  Five  out  of  the  eight 
data  points  of  aircraft  2  were  outside  the  ±1  PR  range. 
Therefore  addition  of  a  constant  bias  did  not  result  in  any 
changes  in  the  correlation. 

The  above  analyses  reveal  that  a  trend  does  exist  uetween 
the  pilot  rating  and  the  pitch  phase  angle  data  obtained  by 
generating  the  open-loop  frequency  response.  However,  the 
pilot  may  not  desire  to  operate  all  the  systems  at  3  rad/sec. 
Therefore,  fixing  the  cross-over  frequency  at  3  rad/sec  may 
not  necessarily  be  an  appropriate  assumption  for  performance 
representation  and  the  workload  evaluation. 

Alternate  Analysis  For  Data  Correlation 

It  seems  reasonable  and  practical  that  the  pilot,  does 
not  close  the  loop  at  3  rad/sec  for  every  case,  but  adjusts 
for  a  suitable  cross-over  frequency  by  adding  some  lead  and 
maintaining  a  desirable  phase  margin.  Pilots  generally  do 
not  object  to  adding  leads  of  the  order  of  0.5  sec  to  1.0 
sec  time  constant.  The  phase  angle  data  of  Table  1  thru  5 
for  the  pitch  attitude  show  that  addition  of  a  pilot  lead  of 
0.5  sec.  to  1.0  sec.  in  general,  would  provide  sufficient 
phase  margin  for  all  aircraft  with  no  control  system 
(C/S  =1)  when  evaluated  at  a  frequency  1  rad/sec  above  the 
aircraft  short  period  natural  frequency.  This  fact  is  also 
shown  by  the  closed  loop  analog  simulation  results,  discussed 
in  the  next  chapter.  However  with  the  addition  of  control 
systems,  the  pilot  tends  to  adjust  the  lead  time  constant 
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and  the  cross-over  frequency  to  obtain  a  desirable  phase  mar¬ 
gin.  As  an  example  to  elaborate  this  point,  let  us  examine 
the  phase  angle  of  all  aircraft  with  no  control  system  and 
with  a  pilot  lead  time  of  either  0.5  or  1.0  sec.  The  phase 
angle  generated  at  the  frequencies  of  interest  and  the  resul¬ 
tant  phase  margins , are  tabulated  in  Table  1. 


Table  1  Phase  Angle  Due  to  Pilot  Lead  at  and  Range  of 


Air¬ 

craft 

No. 

WSP 

— 

Uc  =  USP+1 

*° 
at  u 

c 

For 

tl=0’5 

^  Lead 

at  u) 

c 

For 

TL=0.75 

‘f’  Lead 

at  a) 

c 

Range  of 
*M  f0r 

0.5  <_  T 

<  1.0 

M 

2 

-196 

+45 

+63 

29  -  67 

^9 

^B 

3.3 

-206 

+59 

+73 

34  -  48 

^9 

3 

-220 

+56 

+72 

16  -  32 

n 

1.99 

3 

-188 

+56 

+72 

48  -  54 

IB 

3.9 

5 

-227 

+63 

+76 

26  -  29 

The  data  shown  in  Table  1  indicate  that  a  sufficient  range 
of  phase  margin  is  available  and  hence  it  may  be  reasonable  to 
choose  uc  a  +  l  .  The  effect  of  adding  control  systems 

to  the  aircraft  in  terms  of  phase  angle  and  slope  at  wc  can 

then  be  measured.  The  overall  concept  of  fixing  the  perfor¬ 
mance  and  correlating  pilot  ratings  with  workload  and  sensi¬ 
tivity  remains  the  same,  except  that  wc  for  each  aircraft  is 


chosen  to  be  (ugp  +  ^  5  insteac*  °f  a  fixed  3  rad/sec.  The 

method  of  least  squares  is  used  to  curve  fit  the  data  for 
correlating  the  pilot  ratings  with  the  respective  phase  angles 
and  slopes  at  <*)c. 

The  analysis  is  not  repeated  for  aircraft  3  and  4  as  it 
has  already  been  demonstrated  earlier  since  for  these  air¬ 
craft  (wgp  +1)  is  3  rad/sec.  Therefore,  aircraft  1  and  2 

are  analyzed  using  2  and  3.3  rad/sec  as  u>c  respectively,  to 

test  the  validity  of  the  choice  of  «c .  Only  the  lag  control 

systems  were  considered  to  draw  a  comparison  between  these 
results  and  the  earlier  results.  The  deviation  in  PR  from  the 
best  ratings  were  correlated.  The  results  thus  obtained  are 
tabulated  in  Table  E15. 

Comparison  of  the  Two  Methods 

The  pilot  rating  expressions  obtained  by  both  methods 
showed  a  similar  trend  between  the  ratings  and  the  phase 
angle  and  slope  at  the  chosen  cross-over  frequency.  Unfor¬ 
tunately,  there  was  no  significant  improvement  when  the  second 
method  was  used.  However  it  only  seems  appropriate  and  reason¬ 
able  to  select  different  cross-over  frequencies  for  each  air¬ 
craft.  The  same  concept  is  used  for  analyzing  the  flight  path 
angle  data. 

Analysis  of  Flight  Path  Angle, n  Data 

The  data  recorded  in  Tables  E6  thru  E10  for  the  phase 
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angle  and  slope  of  phase  curve  of  ' y '  exhibit  different  char¬ 
acteristics  than  the  pitch  angle  data.  Since  the  transfer 

e 

function  constitutes  a  non-minimum  phase  system,  the  system 

0 

tends  to  be  unstable  much  faster  than  the  minimum  phase  — 

e 

transfer  function.  From  the  data  it  is  seen  that  if  the  pilot 
closes  the  y-loop  at  a  frequency  1  to  1.5  octave  below  the  w^p, 

he  is  able  to  operate  the  aircraft  with  no  control  system 
dynamics  with  a  nominal  0.5  to  1.0  sec  lead  time  constant. 
Therefore,  instead  of  analyzing  the  data  at  a  fixed  the 

phase  angle  and  slope  of  the  phase  curve  at  a  frequency  1 
octave  below  Wgp  was  analyzed  to  match  a  pilot  rating  expres¬ 
sion  of  the  form 


PR  -  “1*  +  “2  I! 


+  a 


6  <{> 
3  6  to 


(15) 


where  a^,  and  are  constants  evaluated  for  each  air¬ 
craft  individually  using  the  method  of  least  squares  explained 


i  •  6  4> 

earlier.  — 

6  w 


is  the  slope  of  the  phase  curve  below  u>c  and 


and  -7!  is  the  slope  of  the  phase  curve  above  to  . 

4>U)  0  r  r  c 


The  values  of  a^>a2  anc*  “3  at1<^  the  estimated  and  actual 


PR  are  tabulated  in  Table  E16.  For  aircraft  1  and  5  the  two 
slopes  are  averaged,  since  the  variation  is  small.  A  plot  of 
the  estimated  PR  versus  the  actual  PR  is  shown  in  Fig.  5. 


ESTIMATED  PR 


Fig.  5.  Estimated  Verses  Actual  PR  Using  V  Phase  Angle 
Data  And  Equation  15 
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The  lead  control  system  configurations  were  also  included  for 
this  analysis.  Out  of  42  data  points  13  data  points  are  out¬ 
side  the  ±1PR  boundary.  8  out  of  the  13  points  represent  air¬ 
craft  1,  and  3  points  represent  aircraft  2,  and  2  points  re¬ 
present  aircraft  4.  If  the  data  for  aircraft  1  is  removed 
from  the  figure,  only  5  points  out  of  32  are  outside  the  ±1PR 
boundary,  which  is  a  much  better  correlation  compared  to  the 
pitch  angle  data. 

To  observe  the  effect  of  adding  control  systems  on  pilot 
rating  an  alternate  expression  was  developed  as  shown: 

APR  =  a. A$  +  a0  A  (16) 

1  Zoo) 

where  APR  is  the  deviation  in  pilot  rating  from  the  best 
rating  in  the  data.  A 4)  and  A  ~  are  the  deviations  in  the 
corresponding  phase  angle  and  slope.  Five  expressions  were 
developed,  one  for  each  aircraft.  The  values  of  and 

and  the  estimated  and  actual  APR  are  tabulated  in  Table  E17. 
Fig.  6  shows  a  plot  of  estimated  versus  actual  APR.  It  is 
seen  that  14  out  of  37  points  lie  outside  the  ±1aPR  boundary. 
Out  of  the  14  points,  5  represent  aircraft  1,  5  represent  air¬ 
craft  2,  2  points  represent  aircraft  4  and  2  points  represent 
aircraft  3  and  5. 

Summary  of  Open-Loop  Frequency  Response  Analysis 

The  results  obtained  from  the  open-loop  frequency 
response  analysis  can  be  summarized  as  follows: 
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Fig. 


Estimated  Verses  Actual  APR  Using  ' y '  Phase  Angle 
Data  And  Equation  16 


1)  For  a  particular  cross-over  frequency  (fixed 
performance)  the  phase  angle  (pilot  workload)  and  the 
slope  of  the  phase  curve  (system  sensitivity)  correlated 
reasonably  well  with  the  pilot  rating.  This  was  especi¬ 
ally  true  for  the  flight  path  angle  data  (Figure  5)>  In¬ 
clusion  of  the  system  performance  as  a  variable  in  the 
pilot  rating  expression  may  provide  better  correlation. 

An  approach  for  accomplishing  this  is  described  in  the 
recommendation  section  of  this  report. 

2)  The  deviation  in  the  pilot  rating  with  respect  to 
the  best  rating  correlated  well  with  the  corresponding 
deviation  in  the  phase  angle  and  slope.  This  was  true 
for  both  the  pitch  angle  (Figure  4)  and  the  flight  path 
angle  (Figure  6)  data. 

3)  The  pilot  can  obtain  good  phase  margin  for  aircraft 
with  no  control  system  dynamics  if, 

a)  He  closes  the  pitch  loop  at  uc  approximately 

equal  to  (“gp  +  1)  rad/sec. 

b)  He  closes  the  y-loop  at  w  approximately  one 
octave  below  w^p. 

4)  Aircraft  with  control  systems  introducing  significant 
lags  around  Wgp  will  have  poor  ratings. 

5)  The  degree  of  correlation  between  the  phase  angle 
data  and  the  pilot  ratings  indicate  that  the  pitch  and 
flight  path  angles  seem  to  be  important  tracking  para¬ 
meters  for  pilots  during  landing  task. 


IV  Analog  Computer  Simulation 

The  Analog  computer  is  a  tool  which  helps  in  the  simu¬ 
lation  and  study  of  a  physical  system.  It  can  be  used  for 
the  analysis,  design  and  testing  of  open-loop  and  closed-loop 
control  systems  in  the  time  domain  for  standard  inputs. 

Analog  simulation  studies  are  very  handy  for  parameter  ad¬ 
justments  in  a  control  system  to  achieve  a  desired  objective, 
vis-a-vis,  performance  and  time  response. 

This  chapter  discusses  the  simulation  of  the  pilot- 
aircraft  closed-loop  system  on  the  analog  computer.  This 
study  is  undertaken  to  find  an  optimum  pilot  lead  time  which 
yields  the  best  system  performance  for  the  aircraft-control 
system  configurations  simulated.  The  performance  thus  meas¬ 
ured  is  correlated  with  PR.  The  objective  is  two  fold: 

1)  To  see  if  a  correlation  between  performance  and  PR  exists; 

2)  To  see  if  the  pilot  could  improve  the  system  performance 
by  choosing  different  forms  of  pitch  angle  and  flight  path 
angle  feedback.  The  "Lead  only"  form  of  McRuer  model  is 
used  for  the  pilot  as  explained  in  Chapter  II  with  a  reaction 
time  delay  of  0.3  sec.  The  ISE  (integral  of  squared  error) 

is  used  as  the  performance  criterion.  The  pilot  gain,  K^,  and 

lead,  T^,  are  adjusted  to  minimize  ISE  when  the  pilot  performs 
a  particular  tracking  task  for  a  step  command.  Using,  the 
short,  period  approximation  dynamics  of  A/C  1,  first  the  effects 
of  adding  control  systems  on  the  performance  are  studied. 
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The  tracking  task  is  analyzed  in  three  phases:  1)  Inner-loop 
pitch  attitude  (o)  tracking  for  a  step  pitch  command;  2)  outer- 
loop  flight  path  angle  (y)  tracking  for  a  step  y-command;  and 
3)  combined  pitch  and  flight  path  angle  tracking  (inner  and 
outer  loop)  to  minimize  the  error  in  the  flight  path  angle 
for  a  step  y-command.  The  results  thus  obtained  are  described. 
These  three  phases  are  repeated  for  the  remaining  four  aircraft 
using  the  baseline  control  system  (control  system  =  1)  to  ob¬ 
serve  the  performance  variations  between  the  five  aircraft 
configurations.  Since  the  natural  frequency  of  the  elevator 
servo  (75  rad/sec)  and  the  feel  system  (26  rad/sec)  dynamics 
are  very  high,  their  effect  at  the  frequencies  of  interest  are 
assumed  to  be  negligible  and  hence  are  not  included  in  the 
simulation. 

Thus,  the  objective  of  the  analog  simulation  study  is  to 
find  optimum  pilot  parameters  and  which  would  minimize 

the  ISE  for  a  particular  tracking  task.  The  description  of 
the  analog  computer  circuits  used  during  the  simulation  may 
be  found  in  Appendix  D.  References  13,14,  and  15,  were  very 
useful  in  programming  the  simulation. 

Pitch  Attitude  Tracking  Task 

The  block  diagram  used  for  simulation  is  shown  in  Fig.  7. 
Aircraft  No.  1  with  no  control  system  dynamics  was  first 
analyzed.  The  control  systems  which  were  flown  with  air¬ 
craft  1  were  then  simulated  and  the  analysis  was  repeated  for 
each  case. 
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Figure  7  Block  Diagram  of  Pitch  Tracking  Task 

The  proceedure  used  for  the  analysis  is  as  follows.  The 

error  signal  e0,  was  fed  to  a  quarter-square  multiplier 

o 

circuit  to  generate  e0,  which  was  then  integrated  to  gener¬ 
ate  ISE.  The  pilot  lead  time  TL  was  varied  from  0  to  1.25 
sec  in  steps  of  0.25  sec.  For  each  value  of  T^,  the  pilot 

gain  K  was  varied  and  the  value  of  ISE  was  recorded  for  a 
P 

step  pitch  command  (6C)- 

The  data  thus  generated  is  shown  in  Fig.  8  thru  16. 

While  pilot  lead  times  greater  than  0.5  or  0.75  sec,  produce 
a  slight  improvement  in  the  system  performance,  they  do  so 
at  the  expense  of  greater  pilot  workload.  For  a  1.0  sec.  and 
1.25  sec.  lead  time  the  relative  difference  in  the  ISE  value 
is  insignificant.  Besides,  the  curves  for  1.0  and  1.25  sec. 
lead  time  show  much  larger  variations  in  the  value  of  ISE 
about  the  minimum,  thus  indicating  higher  system  sensitivity 
due  to  pilot  gain  variations.  Therefore,  a  value  of  0.5  or 
0.75  sec.  pilot  lead  seems  to  be  the  optimum  choice.  Fig.  14 


41 


shows  the  performance  curve  of  aircraft  1  with  control  system 
2 

dynamics  -  .  Since  this  control  system  introduces 

(S  +  2) 

significant  amount  of  lag  near  the  short  period  natural  fre¬ 
quency  of  the  aircraft  (1  rad/sec),  the  aircraft  dynamic  re¬ 
sponse  is  very  sensitive  to  gain  change.  In  other  words,  for 
good  performance  the  pilot  must  hold  the  gain  steady  to  avoid 
large  variations  in  and  4^.  This  in  effect  amounts  to  an 

additional  workload  in  terms  of  pilot  concentration  and 
hence  a  degraded  pilot  opinion  of  the  task.  This  corresponds 
to  a  very  negative  slope  of  the  phase  curve  around  wc  in  the 

open-loop  frequency  response  analysis  of  the  last  chapter. 

The  value  of  minimum  ISE  and  the  corresponding  pilot  ratings 
are  plotted  in  Fig.  17  and  18  for  a  0.5  ard  0.75  sec.  pilot 
lead,  respectively.  The  rating  degradation  with  poor  per¬ 
formance  for  lag  cases  can  be  observed,  but  the  reverse  trend 
with  the  lead  cases  is  again  a  puzzle  which  could  not  be 
explained . 

Flight  Path  Angle  Tracking  Task 

During  the  flare  maneuver,  the  sink  rate  information  is 
one  of  the  visual  cues  the  pilot  uses  in  performing  the 
landing  task.  Since  flight  path  angle,  y  ,  is  directly  related 
to  the  sink  rate,  it  is  logical  to  analyze  the  pilot  in  the 
•y-loop.  This  was  done  by  first  assuming  no  pitch  angle  feed¬ 
back.  The  block  diagram  used  for  this  simulation  is  shown  in 
Fig.  19.  Short  period  dynamics  of  aircraft  1  were  used  and 
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the  ISE  (  e^)  was  measured  for  a  step  y-command.  The  pilot 

parameters  K  and  were  varied  and  the  procedure  used  for 
PY 

the  analysis  of  pitch  loop  was  repeated,  vis-a-vis,  measuring 

the  value  of  ISE  until  a  minimum  is  reached,  as  K  and  T. 

PY 


were  varied. 


Figure  19  Block  Diagram  of  y  Tracking  Task 

The  plots  of  ISE  versus  the  gain,  K  as  T^  is  varied  is 

PY 

shown  in  Figures  20  thru  28.  From  the  plots  it  is  seen  that 
the  r-loop  performance  shows  gradual  improvement  as  pilot  lead 
time,  T^,  increases.  However,  the  relative  improvement  in 

performance  beyond  T^  =  0.75  sec.  shows  a  converging  trend 

and  indicates  that  additional  pilot  lead  beyond  this  point 
does  not  result  in  proportional  improvement  in  performance. 
Therefore  a  0.75  sec.  lead  seems  to  be  the  optimum  lead  time 
constant  for  the  closed-loop  flight  path  angle  tracking  task. 
The  tendency  of  the  system  to  become  more  and  more  sensitive 
in  the  presence  of  lag  control  systems  is  evident  from  the 


Fig.  22  ,  ISE  Vs.  K  Aircraft  1 
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Fig.  25  ISE  Vs.  K  -  Aircraft  1  C/S  5 

Pr 
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plots  when  exceeds  0.75  sec.  The  addition  of  lag  control 

system  dynamics  also  results  in  a  system  which  is  more  sensi¬ 
tive  to  gain  variations.  In  particular  Fig.  26,  the  perfor¬ 
mance  curve  of  aircraft  1  with  control  system  6,  2/S+2,  depicts 
a  very  sensitive  system.  The  achievable  minimum  ISE  is  also 
considerably  greater  than  the  minimum  ISE  for  the  aircraft 
with  the  other  control  systems.  Similar  results  for  this 
aircraft /control  system  combination  were  observed  in  the 
pitch  tracking  task  and  in  the  open-loop  frequency  response 
analysis  of  Chapter  III.  Figures  29  and  30  show  a  plot  of 
minimum  ISE  versus  the  PR.  The  plots  show  an  increase  in 
rating  as  the  performance  degrades.  The  reverse  trend  for 
the  lead  control  systems  needs  further  investigation. 


Pitch  and  Flight  Path  Angle  Tracking 


The  above  analyses  show  that  both  pitch  and  flight  path 
angles  are  important  parameters  for  the  pilot  during  the  land¬ 
ing  task.  The  next  logical  step  is  to  study  the  effect  of 
feeding  Lack  both  pitch  and  flight  path  angle  as  a  multi¬ 
output  system  for  a  single  input  step  y-command.  Figure  31 
shows  a  block  diagram  of  the  system  used  for  simulation. 


Fig.  31  Both  0  and  y  Tracking  Task  For  A  Step  T-Command 
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The  analysis  was  restricted  to  only  aircraft  1  and  the 
lag  control  systems  flown  with  aircraft  1.  The  short  period 
configuration  was  used  for  the  aircraft  simulation.  The 
'lead  only'  form  of  pilot  dynamics  was  used  for  the  pilot  in 
the  pitch  loop.  Two  values  of  T^,  0.5  and  0.75  seconds,  were 


used  and  the  corresponding  gains  K 


which  minimized  the  ISE 


in  the  pitch  tracking  task  were  used.  With  the  inner  loop 

parameters  fixed,  the  closed-loop  pitch  transfer  function 

forms  a  type  zero  system.  If  only  a  pilot  gain,  K  ,  is  used 

r 

in  the  outer  loop  there  will  always  be  a  steady-state  error, 
e^ ,  for  a  step  y-command  input.  Therefore,  an  integral  plus 
proportional  control  system  was  used  as  shown  in  Figure  31. 

It  is  assumed  that  the  integral  plus  proportional  control 
represents  the  pilot  in  the  outer  loop.  A  pilot  reaction 
time  delay  was  not  included  in  the  outer  loop. 

The  transfer  function  used  for  integral  and  proportional 
control  is  of  the  form; 


G(s)  =1+| 


(17) 


Determining  the  appropriate  value  of  'a'  to  be  used  for  the 
» 

task  was  the  next  logical  step.  Therefore,  K  was  varied 

PY 

2 

to  measure  ISE  /e^  for  a  step  y-command  while  the  value  of 

'a'  was  varied  from  0.1  to  0.4  in  steps  of  0.05.  This  pro- 
ceedure  was  used  for  both  a  T^  of  0.5  and  0.75  seconds  for 

pilot  lead  in  the  pitch  loop.  The  data  thus  generated  is 
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shown  in  Tables  2  and  3.  Since  the  change  in  the  value  of  the 
ISE  for  'a'  greater  than  0.25  was  not  significant,  a  value  of 
a  =  0.25  was  selected  and  used  for  subsequent  analysis. 


The  lag  control  systems  flown  with  aircraft  1  were  sim¬ 
ulated  to  study  the  effect  of  performance  degradation  on 
pilot  rating.  The  values  for  the  variables  Kp  ,  a  ,  and  T^ 


were  discussed  above.  The  outer  loop  pilot  gain  K  was  var- 

PY 

2 

ied  and  the  ISE  /e^  was  recorded  for  a  step  Y-command. 

2 

Figures  32  and  33  show  a  plot  of  ISE  /e  versus  the  K  .  The 

Y  PY 


effect  of  performance  degradation  on  pilot  rating  can  be  clear¬ 
ly  seen  from  these  plots.  Figures  34  and  35  show  the  plot  of 
o 

minimum  ISE  /ey  versus  the  PR  for  aircraft  1  without  the  con¬ 
trol  system  and  with  the  lag  control  systems.  These  figures 
also  indicate  a  relationship  between  PR  and  performance  for 
a  fixed  pilot  workload. 


Comparison  of  Flight  Path  Angle  Tracking  Task  With  and  With¬ 
out  the  Inner  Pitch  Loop 

The  flight  path  angle  tracking  without  the  inner  pitch 

2 

loop  yields  a  relatively  higher  value  of  ISE  /e  as  com¬ 
pared  to  the  tracking  with  the  pitch  loop.  However,  the  pilot 
model  becomes  more  complicated  for  the  two  loop  closure  task, 
in  that,  in  addition  to  the  inner  loop  lead  time  and  gain 
which  is  common  to  both  the  tasks,  an  outer  loop  gain  and  in¬ 
tegral  plus  proportional  control  is  required  from  the  pilot. 
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0.15  a  =  0.2  |  a  =  0.25  a  =  0.3  a  =  0.35  a  =  0.6 
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Table  3  K  Verses  ISE  For  Tt  =  0.75  -  Different  Values  Of  "a’ 
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Fig.  34  PR  Vs.  ISE)j>jin  For  =  0.5  Aircraft  1  And  Lag  C/S 
*6’  And  'v 1 ' TrarVinp  TncW 


The  task  of  integrating  an  error  essentially  means  remember¬ 
ing  the  past  errors  and  thus  represents  an  increase  in  the 
pilot  workload.  Whether  the  pilot  would  like  this  kind  of 
additional  workload  or  if  the  additional  workload  is  worth 
the  improvement  in  the  performance  needs  to  be  determined. 
Also,  for  a  multi-output  tracking  task,  such  as  the  one  de¬ 
scribed  here,  how  the  pilot  would  time  share  his  control  of 
the  two  outputs  e  and  y  ,  needs  to  be  investigated.  It  may 
be  possible  that  during  the  glide  slope  the  pilot  is  flying 
the  aircraft  by  tracking  only  the  pitch  angle,  e  ,  up  to  a 
certain  altitude  where  the  rate  of  sink,  h  becomes  important 
for  him.  This  altitude  may  be  about  100  feet,  where  the 
pilot  may  begin  tracking  both  pitch  angle  and  flight  path 
angle  until  the  flare  maneuver.  During  flare,  especially 
the  last  50  feet,  the  pilot  may  switch  over  to  tracking  only 
the  altitude,  assuming  that  he  has  done  his  best  to  stabilize 
the  o-loop.  A  plot  of  ISE  versus  Kp  is  shown  in  Figure  36 

for  a  pitch  loop  T^  of  0.5  and  0.75  seconds  for  the  three 
tracking  tasks  using  aircraft  1  without  the  control  systems. 

The  above  hypothesis,  i.e-,  during  approach  the  pilot  is 
essentially  tracking  the  pitch  attitude  and  during  flare  and 
landing  phase  the  pilot  is  tracking  the  flight  path  angle 
while  holding  a  fixed  optimum  pitch  loop  performance  may  be 
verified  by  the  following  proceedure.  The  pilot,  in  a  large 
number  of  flight  configurations  flown,  evaluated  the  approach 
task,  and  the  overall  landing  task  separately  and  provided 
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y  -  Loop  Closed 
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ISE  Vs.  K  /K  Aircraft  1  C/S 
”y  "0 


two  different  ratings.  An  optimum  pilot  workload  and  perfor¬ 
mance  evaluation  can  be  performed  on  analog  computer.  The 
proceedure  described  in  this  chapter  and  those  flight  con¬ 
figurations  for  which  both  ratings  were  provided  may  be  used; 
first  as  a  pitch  tracking- task,  and  then  as  a  y-loop  tracking 
task  with  the  pitch  loop  optimized  as  per  the  first  analysis. 
The  results  thus  obtained  can  be  correlated  with  the  corres¬ 
ponding  approach  and  overall  pilot  ratings  to  validate  the 
hypotheses.  This  analysis  was  not  carried  out  due  to  lack  of 
time.  The  data  generated  in  this  chapter  is  insufficient  to 
draw  any  positive  conclusions  about  this  hypothesis. 

Pitch  Attitude  Tracking  For  Different  Aircraft 

The  previous  analysis  is  restricted  to  only  Aircraft  1. 
In  order  to  extend  the  analysis  to  aircraft  configurations 
2,3,4,  and  5,  it  is  assumed  that  the  effects  of  adding  con¬ 
trol  systems  will  have  similar  trends  as  observed  for  air¬ 
craft  1.  Therefore  only  the  aircraft  short  period  dynamics 
are  simulated  and  analyzed  without  adding  the  control  systems. 
The  analysis  is  carried  out  in  two  phases;  1)  closing  the 
pitch  loop  alone  for  a  step  ©-command.  2)  closing  the 
y-loop  without  closing  the  pitch  loop  for  a  step  y-command . 

A  pilot  lead  time  of  0.5  sec.  and  0.75  sec.  was  selected  and 
the  pilot  gain  Kp  was  varied  for  each  T^  to  measure  ISE  for 

a  step  0-command.  The  choice  of  0.5  sec.  and  0.75  sec.  lead 
time  is  based  on  the  previous  analysis  of  aircraft  1.  The 
values  of  K  and  ISE  thus  measured  for  each  aircraft  are 


plotted.  Figure  37  shows  a  plot  of  ISE  versus  the  pilot 
gain  for  the  five  aircraft  configurations  with  a  lead  time 
of  0.5  sec.  The  figure  shows  that  aircraft  1,3  and  5  are 
relatively  more  sensitive  to  gain  variations  than  aircraft 
2  and  4.  The  values  of  and  phase  margin  (<t>^)  were  de¬ 
termined  from  the  Bode  plot  for  the  gain,  K  ,  corresponding 

Pe 

to  the  minimum  value  of  ISE.  Digital  computer  program  TOTAL 
was  used  for  this  purpose.  Table  4  shows  the  values  of  “ 

and  <t» for  each  aircraft  configuration. 

The  values  of  minimum  ISE  for  each  aircraft  configura¬ 
tion  are  plotted  against  PR  (Figure  38).  While  the  plot 
does  not  show  a  direct  relationship  between  PR  and  perfor¬ 
mance  (ISE),  a  rating  expression  of  the  form  shown  in  Eqn.  18 
indicates  that  a  linear  combination  of  —  and  -r~  has  a 


c. 


M 


relationship  with  PR. 


P.R.  =  a.  -i- 
1  (o 

c 


+  a. 


M 


(18) 


The  values  of  and  <»2  were  determined  using  the  method  of 

least  squares.  the  plot  of  actual  versus  estimated  PR  using 
Eqn.  18  is  shown  in  Figure  39.  The  above  proceedure  was  re¬ 
peated  for  a  0.75  sec.  lead  time  and  the  data  thus  generated 
is  plotted.  A  plot  of  ISE  versus  is  shown  in  Figure  40. 

From  the  figure  it  can  be  clearly  seen  that  all  aircraft  are 
relatively  more  sensitive  to  gain  variations  as  compared  to 


74 


Fig.  38  PR  Vs.  ISE)Hin  Different  Aircraft  Pitch  Tracking 
Task  T^  =  0.5  sec. 
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Aircraft  No. 

0) 

C 

rad/ 
sec . 

Pilot 

Ratings 

Est . 

Actual 

1 

1.8 

40 

3.38 

4 

2 

3.5 

40 

3.03 

2 

3 

3 

31 

3.98 

4 

4 

3 

60 

2.2 

2 
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5.5 
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4.5 

5 
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Table  4  u<  And  For  Different  Aircraft  Evaluated 
c  M 

At  ISE)^n  And  =  0.5  For  Pitch  Tracking  Task 


Aircraft  No. 

0) 

C 

rad/ 
sec . 

♦m 

Pilot 

Ratings 

Est . 

Actual 

1 

2.1 

45 

3.46 

4 

2 

4 

37 

3.1 

2 

3 

3.3 

30 

3.78 

4 

4 

3.4 

60 

2.38 

2 

5 

6 

20 

4.68 

5 

a1=3.42 

a2=82.27 

Table  5  to  And  <t>w  For  Different  Aircraft  Evaluated  At 
c  M 

ISE)M.  And  T.  =  0.75  For  Pitch  Tracking  Task 


a  0.5  sec.  lead  time.  The  values  of  wc  and  4>M  obtained  from 

the  Bode  plot  for  the  gain  K  ,  corresponding  to  the  minimum 

pe 

ISE  (Table  5),  show  that  with  a  higher  lead  time  the  pilot  is 
able  to  close  the  pitch  loop  at  a  slightly  higher  wc  at  the 

cost  of  increased  sensitivity.  A  plot  of  PR  versus  minimum 
ISE  for  the  five  aircraft  (Figure  41)  for  0.75  sec.  lead  time 
shows  that  four  out  of  five  points  are  within  ±1  PR  limits. 

The  estimated  PR  were  evaluated  using  Eqn.  18.  A  plot  of 
actual  and  estimated  PR  for  a  0.75  lead  time  using  Eqn.  18  is 
shown  in  Figure  42.  From  the  figure  it  can  be  seen  that 
1  1 

—  and  —  have  a  proportional  relationship  with  the  PR. 
u  c 

The  values  of  w  obtained  for  the  five  aircraft  for  a  0.5 
c 

sec.  lead  time  and  for  minimum  ISE  indicate  that  the  pitch 
attitude  task  performance  is  optimum  for  a  value  of  wc  which 

is  approximately  one  rad/sec  higher  than  Wgp.  A  similar  re¬ 
sult  was  obtained  in  Chapter  III.  Another  conclusion  revealed 
from  the  open-loop  frequency  response  analysis  is  the  closure 
of  y— loop  at  a  an  octave  below  Wgp  .  To  confirm  this  result, 

the  closed-loop  flight  path  angle  tracking  task  was  simulated 
as  discussed  below. 

Flight  Path  Angle  Tracking  Task-Different  Aircraft 

The  proceedure  used  for  the  analysis  of  closed-loop 
flight  path  angle  tracking  task  was  the  same  as  for  the  pitch 


^ 


angle  tracking  task.  The  values  of  were  adjusted  to  0.5 

and  0.75  sec.  to  obtain  two  sets  of  readings.  The  value  of 

pilot  gain  K  was  varied  to  obtain  ISE  for  a  step  y-command. 
PY 

The  plot  of  K  versus  ISE  is  shown  in  Figures  43  and  44,  for 

PY 

0.5  sec.  and  0.75  lead  time  respectively. 

Figures  43  and  44  reveal  that  aircraft  1  and  3  are  more 
sensitive  to  gain  variation  than  aircraft  2,4,  and  5  for  the 
Y-tracking  task.  The  increase  in  sensitivity  with  increase 
in  pilot  lead  time  can  also  be  observed.  The  values  of 

and  were  determined  for  Kp  corresponding  to  minimum  ISE 

value  and  tabulated  in  Tables  6  and  7.  The  values  of  «c  thus 
obtained  show  that  for  aircraft  1,2  and  4,  is  very  nearly 
an  octave  below  while  for  aircraft  5,  it  is  1.5  to  2  oc¬ 


taves  below  wcp«  Aircraft  3,  however,  yields  of  2  rad/sec 

JSP 


SP‘  -  - ,  y - -c 

which  is  the  value  of  o)Qp  and  thus  contradicts  this  guideline. 


A  plot  of  estimated  versus  actual  PR  is  shown  in  Figures  45 
and  46  using  Eqn.  18  for  T^  of  0.5  and  0.75  sec.  The  plots 


show  that  for  aircraft  1,2,3  and  4  a  correlation  exists,  but 
the  estimated  rating  of  aircraft  5  needs  further  investiga¬ 
tion.  It  appears  that  PR  is  a  function  of  only  — -  within 


M 


certain  limits.  According  to  this  relationship,  higher  values 
of  will  reduce  the  rating,  but  we  know  that  the  system  per¬ 
formance  for  greater  than  80°  will  be  sluggish  and  the 
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•’lg.  43  K  Vs.  ISE  y-  Tracking  Task  -  Different  Aircraft  For 


Tracking  Task  Different  Aircraft  For  T,  ■  0.75 


Aircraft  No. 


0) 

C 

i  0 

*  M 

Pilot 

Ratings 

rad/ 
sec . 

Est . 

Actual 

0.6 

43 

4.96 

4 

1.0 

66 

2.04 

2 

1.9 

40 

3.49 

4 

0.9 

45 

3.01 

2 

1.0 

80 

2.94 

5 

a*=2.6  a0=26.71 


Table  6  u>c  and  Evaluated  at  ISE)^in  y-Tracking 

Task  Different  Aircraft  With  T^  =  0.5 


(A) 

c 

*  M 

Pilot 

Ratings 

rad  / 
sec . 

Est . 

Actual 

0.7 

43 

4.24 

4 

1.3 

66 

2.6 

2 

2.0 

42 

3.34 

4 

1.1 

47 

3.46 

2 

1.6 

80 

2.1 

5 

a, =1.04  a„=118.53 


Table  7  a>c  And  4> ^  Evaluated  at  ISE)^n  y-Tracking 
Task  Different  Aircraft  With  Tt  =  0.75 


ratings  are  likely  to  increase.  Therfore ,  a  boundary  of  wc 
and  <t>^  needs  to  be  determined  for  this  analysis. 

Summary  of  Results  of  Closed-Loop  Analog  Simulation 

The  analog  simulation  studies  described  in  this  chapter 
were  based  on  finding  the  effects  of  control  systems  on  the 
performance,  while  adjusting  the  pilot  workload  (T^)  for  an 
optimum  value  and  then  holding  it  constant.  Different  air¬ 
craft  were  then  analyzed  on  the  same  lines.  The  analysis  con¬ 
sisted  of  both  pitch  and  flight  path  angle  tracking  tasks. 

The  following  are  the  results  of  the  study: 

1)  A  lead  time  of  0.5  to  0.75  seconds  seems  to  be  the 
optimum  pilot  lead  for  most  of  the  aircraft-control 
system  combinations. 

2)  Control  systems  generating  significant  lag  at  the 
short  period  natural  frequency  tend  to  make  the  system 
very  sensitive  to  gain  changes. 

3)  Performance,  measured  in  terms  of  ISE,  has  a  direct 
correlation  with  PR  for  aircraft  with  lag  control  systems. 

4)  The  pilot  ratings  for  different  aircraft  with  no 
control  system  dynamics  added,  shows  a  correlation  with 

performance  measured  both  in  terms  of  ISE  and  and 

—  .  However,  the  range  of  u  and  <j>M  in  which  this  re- 

4>m  cm 

suit  is  valid  needs  to  be  determined. 


5)  The  pilot,  in  general,  would  achieve  the  best  perfor¬ 
mance  if  he  closes  the  pitch  loop  at  a  frequency  «c  equal 

to  Ugp  +  1  rad/sec  for  the  pitch  tracking  task  alone 

and  uc  equal  to  1  to  1.5  octave  below  Wgp  for  flight 

path  angle  tracking  task  alone  for  all  aircraft  without 
the  addition  of  control  systems. 

6)  Closing  both  the  pitch  and  flight  path  angle  loops 
for  a  y-command ,  improves  the  performance  but  at  the  cost 
of  increased  pilot  compensation.  However,  the  relative 
merits  and  demerits  of  closing  only  the  y-loop  versus 
closing  the  pitch  and  y-loop  need  further  study. 

Some  of  the  conclusions  drawn  from  this  study  do  support 
the  earlier  analysis  of  Chapter  III.  Both  chapters  concen¬ 
trate  on  the  effects  of  performance  and  work  load  on  the  pilot 
ratings.  However,  no  attempt  was  made  in  either  of  the 
analyses  which  would  help  in  understanding  the  pilot  comments 
associated  with  each  aircraft  control  system  flown  (Ref  1). 

The  only  observation  made  during  the  closed-loop  pitch  tracking 
task  is  the  noticeable  initial  delay  in  the  closed-loop  time 
histories.  Some  of  the  representative  time  histories  are 
shown  in  Figures  47  and  48.  The  pilots'  comments  indicate 
that  the  initial  pitch  attitude  response,  oscillatory  nature 
of  aircraft  pitch  response,  and  the  final  (steady  state) 
response  seem  to  be  some  of  the  aircraft  behavior  characteris¬ 
tics  that  the  pilot  is  paying  attention  to  in  the  execution  of 
the  landing  task.  Therefore,  an  open-loop  impulse  response 


e(T) 


Fig.  47  Closed-Loop  Pitch  Response  A/C  1  and  Lead  C/S  For 


C/S  6 


study  is  carried  out  using  aircraft  short  period  pitch 
dynamics  to  correlate  some  of  the  pilot  comments  with  the 
aircraft  pitch  attitude.  The  next  chapter  describes  the  open- 
loop  impulse  response  study. 
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This  chapter  presents  the  analysis  of  the  open-loop 
aircraft  pitch  attitude  response  to  an  impulse  input  and  the 
results  thus  obtained.  Pitch  response  to  an  impulse  input  is 
identical  to  the  study  of  pitch  rate  to  a  step  input.  This 
portion  of  the  analysis  resulted  from  the  comments  provided  by 
the  pilot  after  each  flight.  A  detailed  summary  of  the  com¬ 
ments  can  be  found  in  Ref  2.  From  the  comments  it  appears  that 
the  pitch  attitude  response  was  one  of  the  concerns  of  the 
pilot.  In  particular,  the  initial  pitch  attitude  behavior  seems 
to  have  a  significant  impact  on  the  pilot's  opinion.  It  is  em¬ 
phasized,  however,  that  the  initial  pitch  attitude  will  also 
directly  affect  the  flight  path  angle  and  hence  the  altitude 
and  sink  rate. 

Figure  49  shows  the  block  diagram  of  the  system  used  for 
the  analysis.  For  the  same  reasons  as  discussed  in  Chapter  II, 
the  phugoid  mode  was  ignored  and  the  short  period  aircraft 
equations  were  used. 


Figure  49  Block  Diagram  of  Open  Loop  Aircraft  C/S  Configuration 


The  feel  and  actuator  dynamics  were  included  for  this  analysis 
while  the  pilot  block  was  not  included  since  the  study  was  con¬ 
cerned  with  the  behavior  of  the  system  to  pilot  inputs.  Digi¬ 
tal  computer  program  TOTAL  was  used  to  both  generate  and  anal¬ 
yze  the  impulse  response  plots.  From  the  plots  the  initial 
pitch  response  was  studied  in  terms  of  time  taken  for  the  sy- 
sytem  to  reach  ten  percent  of  the  final  value.  Also  an  attempt 
is  made  to  explain  some  of  the  pilot  comments  based  on  the  im¬ 
pulse  response  plots. 

The  procedure  used  for  the  pitch  response  analysis  was  to 
plot  the  pitch  attitude  to  an  impulse  input  for  all  the  air¬ 
craft-control  system  configurations  flown.  The  selected  plots 
thus  generated  are  shown  in  Appendix  F.  The  time  to  reach  ten 
percent  of  final  value  (T^q)  was  recorded  for  all  the  configur¬ 
ations.  Figures  50  and  51  show  a  plot  of  PR  versus  T^q.  The 

plots  clearly  indicate  that  the  pilots  always  liked  the  basic 
aircraft  configuration  (control  system  =  1)  best.  Addition  of 
either  lead  or  lag  control  system  degraded  his  opinion.  This 
trend  indicates  that  the  pilot  anticipates  a  particular  kind 
of  initial  pitch  response.  Any  deviation  from  this  anticipated 
behavior  surprises  the  pilot  causing  him  to  adjust  to  the 
actual  response  he  has  observed.  This  degree  of  adjustment  on 
his  part  is  reflected  in  his  opinion  and  comments.  Addition 
of  a  lead  control  system,  for  instance,  makes  the  initial  re¬ 
sponse  faster  and  more  abrupt  than  his  expectation.  Addition 
of  lag  control  system  introduces  significant  amount  of  initial 


delay  before  the  system  begins  to  respond  to  the  impulse  input. 
The  Pilot  Comments 

After  each  flight  test  performed  in  Ref  2,  the  pilot  was 
asked  to  comment  on  the  initial  and  the  steady-state  pitch 
attitude  response,  especially  the  initial  delay.  Many  of  the 
comments  made  by  the  pilot  can  be  clearly  observed  in  the  open- 
loop  impulse  response  time  histories.  Some  of  the  comments  and 
the  associated  impulse  responses  are  described  in  the  subse¬ 
quent  paragraphs. 

Aircraft  No.  1 

The  pilot  comments  for  aircraft  1  with  lead  control 
systems  indicate  that  the  "initial  pitch  response  was  reason¬ 
able,  but  the  final  response  tended  to  drift".  The  correspon- 
ing  impulse  response  plots  indicate  that  the  duplicating  time 
(Tp)  is  about  one  second  but  the  settling  time  (Tg)  is  about 
A. 8  seconds  (Figures  F1,F2,F3).  The  pilot  has  the  ability 
to  reduce  the  gain  and  adjust  the  initial  fast  response  to  an 
acceptable  initial  response.  But  the  reduction  in  gain  directly 
affects  the  final  (steady-state)  response.  The  pilot  commented 
about  this  final  response  as  "drifting  tendency".  Therefore, 
when  lead  control  systems  are  added,  if  the  pilot  adjusts  the 
gain  for  a  better  initial  response,  he  notices  a  drifting  final 
response,  he  observes  an  abrupt  and  fast  initial  response.  The 
pilot,  therefore,  has  to  work  for  a  balance  between  the  initial 
and  final  acceptable  response,  which  adds  to  his  workload  and 


degrades  his  rating.  This  tendency  is  very  clearly  observed 
in  Figure  47  (Chapter  IV).  The  lead  control  systems  provid¬ 
ing  faster  initial  response  exhibit  a  reverse  trend  in  the 
final  pitch  response.  This  was  the  only  observation  made 
during  the  study  which  explains  the  degradation  in  pilot 
ratings  when  lead  control  systems  were  added.  The  addition  of 
lag  control  systems  produces  an  initial  delay  in  the  response 
ranging  from  0.1  to  0.15  sec  (Figures  F5,F6,F7).  The  pilot 
complained  about  the  "initial  slow  response  and  unpredictable 
final  or  steady-state  behavior."  As  a  result  of  the  initial 
lag  and  unpredictable  final  behavior,  the  pilot  also  com¬ 
plained  about  the  "tendency  to  PIO"  (Pilot  Induced  Oscillations). 
An  initial  delay,  thus,  causes  poor  predictability  and  PIO. 

Aircraft  No.  2 

The  pilot  comments  for  the  base  line  aircraft  (c/s  =  1) 
and  aircraft  with  control  system  No. 3  indicate  that  the  "pilot 
en  joyed  flying"  these  aircraft  configurations.  The  impulse 
response  plots  show  a  very  good  response  (Figures  F9,F10). 

The  response  settling  time  is  about  3.2  sec.  However,  a  large 
overshoot  exists  which  is  being  ignored  by  the  pilot.  The 
pilot  comments  concerning  control  system  No.l  indicate  an 
"abrupt  initial  response"  and  "lack  of  predictability."  The 
impulse  response  shows  that  the  peak  value  is  reached  in  only 
0.4  sec,  while  it  takes  3.9  sec.  for  the  response  to  settle 
down  (Figure  F8).  This  lack  of  predictability  due  to  addition 
of  lead  control  system  has  been  explained  earlier.  The 
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addition  of  lag  control  systems  clearly  shows  an  initial  delay 
(Figures  F11,F12).  The  pilot  comments  reflect  the  "initial 
delayed  pitch  response"  and  the  associated  "PIO  tendency." 

Aircraft  No.  3 

The  short  period  damping  ratio  of  aircraft  No. 3  was  0.25. 
Since  it  was  an  underdamped  system  the  impulse  responses  were 
slowlv  decaying  oscillations.  The  pilot  in  his  comments 
clearly  indicated  this  dynamic  behavior  by  saying  "oscilla¬ 
tory  pitch  response"  or  "looseness  in  pitch."  One  of  the 
pilots,  for  some  reason,  liked  flying  this  aircraft  (PR  =  2) 
with  a  lead  control  system  (No. 3)  which  is  surprising 
(Figure  FI 3).  However,  this  was  one  odd  case,  as  another 
pilot  rated  the  same  aircraft  as  5.  All  flights  with  lag  con¬ 
trol  systems  were  rated  7  and  above.  For  all  the  lag  cases 
the  pilot  complained  about  the  initial  lag,  poor  predictabil¬ 
ity,  and  PIO  tendency.  After  flying  the  aircraft  with  control 
system  5  the  pilot  even  commented  about  his  inability  to  con¬ 
trol  quickly  and  precisely.  The  corresponding  pitch  response 
(Figure  F15)  shows  a  large  peak  time  (1.1  sec.)  and  settling 
time  (8.87  sec.).  He  also  noted  the  "aircraft  bouncing"  in 
turbulance.  This  was  obviously  dv  a  to  poor  damping  ratio  (?gp)* 

Aircraft  No.  U 

The  damping  ratio  (Sgp)  for  this  aircraft  was  1.06,  an 

overdamped  aircraft.  The  pilot  liked  flying  this  aircraft 
without  any  control  systems  added  (PR  =  2)  .  His  comments 
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indicate  good  initial  response,  no  PIO  tendency,  and  no  problem 
with  predictability.  The  corresponding  impulse  response,  shows 
a  smooth  pitch  behavior  (Figure  F18).  However,  when  lead  con¬ 
trol  system  (no. 3)  was  added  and  flown,  the  pilot  complained 
about  the  initial  response  as  "fast  and  abrupt,"  which  caused 
minor  problems  in  landing.  However,  his  rating  (PR  =  3) 
indicates  that  while  he  is  concerned  about  the  fast  and  abrupt 
initial  response,  he  does  not  seriously  object  to  it  (Figure  F17). 
When  lag  control  systems  were  added  and  flown,  he  did  comment 
about  delayed  initial  response,  reduced  predictability  and  the 
resulting  PIO  tendency.  In  one  of  the  cases  (c/s  =  12)  his 
comment  about  initial  pitch  response  read,  "Felt  like  it  was 
glued  in  concrete,"  (PR  =9)  .  The  corresponding  impulse 

response  (Figure  F20)  shows  that,  there  is  a  0.15  sec.  delay 
before  the  response  begins. 

Aircraft  No.  5 

The  short  period  natural  frequency  for  this  aircraft  was 
3.9  rad/sec  and  the  damping  ratio  was  0.54.  This  aircraft  was 
not  flown  with  lead  control  systems.  The  pilot  comments  for 
the  base  line  aircraft  (c/s  =  1)  indicate  that  the  aircraft 
was  "overly  sensitive."  The  pilot  noticed  an  initial 
"hobbling  tendency"  in  pitch  and  expressed  uncertainty  in  pre¬ 
dictability  of  the  steady-state  response.  One  of  the  pilots 
observed  "Quick,  abrupt  and  objectionable"  initial  response 
and  PIO  tendency  during  flare  and  touch  down.  The  correspond¬ 
ing  pitch  impulse  response  plot  shows  a  0.06  sec.  rise  time 


and  2.7  sec.  settling  time  (Figure  F22).  The  addition  of  lag 
control  systems  did  not  seem  to  cause  a  noticeable  delay  in 
the  initial  response  except  for  control  system  13,  where  a 
0.15  sec.  delay  is  noticeable  followed  by  a  fast  rising  re¬ 
sponse  (Figure  F25).  This  tendency  is  also  reported  by  the 
pilot  as  "initial  response  delayed,  then  a  fast  response." 

The  pilot  ratings  for  all  configurations  varied  from  5  to  7 
which  indicate  that  it  was  a  difficult  aircraft  to  fly.  A 
high  frequency  PIO  was  observed  by  the  pilot  in  all  flight 
configurations . 

Summary  of  Impulse  Response  Analysis 

The  aircraft  pitch  attitude  response  to  an  impulse  input 
revealed  the  following: 

1)  Pilots  seem  to  like  a  particular  transient  response. 

Any  faster  or  slower  behavior  degrades  his  rating.  The 
exact  nature  of  this  initial  pitch  response  could  not  be 
identified  in  this  study. 

2)  Addition  of  lead  control  systems,  while  providing  a 
faster  initial  response,  cause  a  drift  in  the  final  re¬ 
sponse.  This  pitch  attitude  behavior  increases  the 
pilot  workload  and  hence  degrades  his  opinion. 

3)  Initial  delay,  unpredictable  steady-state  behavior, 
and  oscillatory  nature  of  pitch  response  due  to  low  damp¬ 
ing  ratio  seem  to  be  the  causes  of  PIO. 

4)  An  initial  delay,  on  the  order  of  0.1  to  0.15  sec.  in 
the  pitch  response  will  in  general  yield  a  poor  pilot  rating. 
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5)  The  time  taken  for  the  pitch  response  to  reach 
ten  percent  of  the  final  value  is  a  reasonable  indication 
of  the  initial  delay  and  may  be  used  as  an  estimate  in 
predicting  deviations  in  the  expected  pilot  rating  with 
reference  to  a  known  rating. 
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VI  Summary  and  Conclusions 


The  conclusions,  supported  by  the  analysis  of  data 
generated  by  using  the  aircraft-control  system  configurations 
from  reference  2  and  the  three  approaches  described  in 
Chapters  III, IV  and  V  are  summarized.  The  three  approaches 
did  reveal  that  pilot  ratings  are  certainly  a  function  of  sy¬ 
stem  performance,  sensitivity,  and  pilot  workload.  In  the 
first  approach  the  system  performance  in  terms  of  wc  and  <!>^ 

is  held  constant  and  the  pilot  opinion  is  analyzed  with  re¬ 
spect  to  pilot  workload  and  system  sensitivity  measured  in 
terms  of  phase  angle  at  w  and  slope  of  the  phase  curve  at  wc 

respectively.  The  conclusions  drawn  from  this  study  are: 

1)  For  a  particular  cross-over  frequency  (fixed  perfor¬ 
mance)  the  phase  angle  (pilot  workload)  and  the  slope  of 
the  phase  curve  (system  sensitivity)  correlated  reason¬ 
ably  well  with  the  pilot  rating.  This  was  especially 
true  for  the  flight  path  angle  data  (Figure  5).  Inclu¬ 
sion  of  the  system  performance  as  a  variable  in  the  pilot 
rating  expression  may  provide  better  correlation.  An 
approach  for  accomplishing  this  is  described  in  the  recom¬ 
mendation  section  of  this  report. 

2)  The  deviation  in  the  pilot  rating  with  respect  to 
the  best  rating  correlated  well  with  the  corresponding 
deviation  in  the  phase  angle  and  slope.  This  was  true 
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for  both  the  pitch  angle  (Figure  4)  and  the  flight  path 
angle  (Figure  6)  data. 

3)  The  pilot  can  obtain  good  phase  margin  for  aircraft 
with  no  control  system  dynamics  if, 

a)  He  closes  the  pitch  loop  at  wc  approximately 

equal  to  (w^p  +  1)  rad/sec. 

b)  He  closes  the  y-loop  at  wc  approximately  one 
octave  below  “gp* 

4)  Aircraft  with  control  systems  introducing  significant 
lags  around  w<,p  will  have  poor  ratings. 

5)  The  degree  of  correlation  between  phase  angle  data 
and  the  pilot  ratings  indicate  that  the  pitch  and  flight 
path  angles  seem  to  be  important  tracking  parameters  for 
pilots  during  landing  task. 

In  the  second  approach  optimum  values  of  pilot  lead  time 
T^  were  determined  by  analog  simulation  studies  using  ISE  as 

a  performance  critierion.  The  pilot  workload  (T^)  was  then 

held  constant  at  the  optimum  value,  and  the  pilot  opinion  was 
studied  as  a  function  of  performance  of  the  aircraft-control 
system  configurations  for  step  inputs.  The  closed-loop  system 
was  simulated  and  analyzed  as  a  pitch  tracking,  flight  path 
angle  tracking,  and  both  pitch  and  flight  path  angle  tracking 
tasks.  The  conclusions  drawn  from  the  results  of  the  Simula- 


1)  A  lead  time  of  0.5  to  0.75  seconds  seems  to  be  the 
optimum  pilot  lead  for  most  of  the  aircraft-control  sy¬ 
stem  combinations. 


2)  Control  systems  generating  significant  lag  at  the 
short  period  natural  frequency  tend  to  make  the  system 
very  sensitive  to  gain  changes. 

3)  Performance,  measured  in  terms  of  ISE,  has  a  direct 
correlation  with  PR  for  aircraft  with  lag  control 
systems . 

4)  The  pilot  ratings  for  different  aircraft  with  no 
control  system  dynamics  added,  shows  a  correlation  with 
performance  measured  both  in  terms  of  ISE  and  -  and 

~  .  However,  the  range  of  u>c  and  <j>^  in  which  this 

^  M 

result  is  valid  needs  to  be  determined. 


5)  The  pilot,  in  general,  would  achieve  the  best  perfor¬ 
mance  if  he  closes  the  pitch  loop  at  a  frequency  coc  equal 

to  «gp  +  1  rad/sec  for  the  pitch  tracking  task  alone  and 

u>c  equal  to  1  to  1.5  octave  below  w^p  for  flight  path 

angle  tracking  task  alone  for  all  aircraft  without  the 
addition  of  control  systems. 

6)  Closing  both  the  pitch  and  flight  path  angle  loops 
for  a  Y-command,  improves  the  performance  but  at  the 
cost  of  increased  pilot  compensation.  However,  the 
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relative  merits  and  demerits  of  closing  only  the  y-loop 
versus  closing  the  pitch  and  r-loop  need  further  study. 

In  the  third  approach  the  aircraft  pitch  attitude 
response  to  an  impulse  input  was  studied  using  the  aircraft- 
control  system  combinations  flown.  The  conclusions  drawn 
from  this  study  are: 

1)  Pilots  seem  to  like  a  particular  transient  response. 

Any  faster  or  slower  behavior  degrades  his  rating.  The 
exact  nature  of  this  initial  pitch  response  could  not  be 
identified  in  this  study. 

2)  Addition  of  lead  control  systems,  while  providing  a 
faster  initial  response,  cause  a  drift  in  the  final  re¬ 
sponse.  This  pitch  attitude  behavior  increases  the  pilot 
workload  and  hence  degrades  his  opinion. 

3)  Initial  delay,  unpredictable  steady-state  behavior, 
and  oscillatory  nature  of  pitch  response  due  to  low  damp¬ 
ing  ratio  seem  to  be  the  causes  of  PIO. 

4)  An  initial  delay  on  the  order  of  0.1  to  0.15  sec.  in 
the  pitch  response  will  in  general  yield  a  poor  pilot 
rating. 

5)  The  time  taken  for  the  pitch  response  to  reach  ten 
percent  of  the  final  value  is  a  reasonable  indication 

of  the  initial  delay  and  may  be  vised  as  an  estimate  in  pre¬ 
dicting  deviations  in  the  expected  pilot  rating  with 
reference  to  a  known  rating. 


VII  Recommendations 


The  following  recommendations  are  made  for  further  study. 

1)  The  first  two  approaches  used  for  analysis  assumed  a 
fixed  performance  and  fixed  workload  respectively.  The  results 
did  reveal  that  A)  correlation  exists  between  PR  and,  system 
sensitivity  and  pilot  workload  (open-loop  frequency  response 
anslysis)  B)  correlation  exists  between  PR  and  performance 
(closed-loop  analog  simulation).  An  alternate  approach,  which 
is  essentially  combining  the  two,  is  to  first  simulate  each 
aircraft-control  system  configuration  flown  on  the  analog  com¬ 
puter  to  determine  optimum  values  of  T^  and  K  .  Then  determine 
the  phase  angle  and  slope  at  wc  from  the  open-loop  frequency 

response  analysis.  Thus  for  each  aircraft-control  system  con¬ 
figuration  there  will  be  a  corresponding  set  of  values  for 
T^,  u>c  and  and  H/fiw,  representing  the  pilot  workload, 

system  performance  and  sensitivity  respectively.  A  linear  re¬ 
gression  analysis  may  then  be  carried  out  between  the  pilot 
ratings  and  the  sets  of  values  for  T^,  o>c  and  <t>^,  s<t>/6u>.  As 

an  alternate  the  measured  value  of  ISE)^n  may  be  used  for 

representing  performance  instead  of  o>c  and  <t>^.  The  procedure 

may  be  used  for  pitch  tracking  task,  flight  path  angle  track¬ 
ing  task,  and  both  pitch  and  flight  path  angle  tracking  task 
simulations  as  explained  in  Chapter  IV.  The  results  thus 
obtained  may  provide  a  better  understanding  of  the  overall 
landing  task. 
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2)  The  pilot  while  evaluating  the  landing  task  provided 
two  separate  ratings:  one  for  the  overall  task  and  one  for 
the  approach  task.  The  data  generated  above  may  be  correlated 
with  each  set  of  ratings  to  see  if  pitch  tracking  task  is  re¬ 
presentative  of  the  approach  task  and  flight  path  angle  track¬ 
ing  task  with  or  without  the  pitch  loop  closed  is  representa¬ 
tive  of  the  overall  task  including  flare  and  touchdown. 

3)  The  range  of  and  4^  needs  to  be  determined  to 

measure  performance  in  terms  of  —  and  —  .  An  alternate 

c  M 

performance  measure  in  terms  of  —  and  4>^  may  also  be  inves- 

wc 

tigated . 

4)  The  exact  nature  of  initial  pitch  transient  response 
and  steady-state  character  liked  by  the  pilots  needs  to  be 
investigated . 

5)  The  degradation  in  pilot  ratings  when  lead  control 
systems  are  added  needs  further  study 

6)  While  analyzing  Aircraft  No. 5,  introducing  a  pilot 
lag  time  constant  (Tj)  to  reduce  the  crossover  frequency  may 

be  considered. 
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Appendix  A 


Aircraft-Control  System  Configurations 
And  Transfer  Functions 


This  appendix  contains  the  transfer  functions  of  the 
thirteen  control  systems  used.  The  aircraft-control  system 
configuration  flown  and  the  corresponding  pilot  ratings  are 
also  listed.  The  first  number  in  the  configuration  represents 
the  aircraft  and  the  second  number  indicates  the  control 
svstem. 

The  feel  svstem  and  actuator  dynamics  were  held  constant 
for  all  configurations.  Approximate  transfer  functions  for 
the  feel  svstem  and  the  actuator  dynamics  are: 

Feel  svstem: 

84.5 

— 

S  +  31. 2S  +  676 

Elevator  Actuator: 

5625 _ 

S2  +  105S  +  5625 

A  description  of  these  two  sub-systems  is  given  in 
Ref.  2. 
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Table : A1  Control  System  Transfer  Functions 


Control 

System 

No: 


Transfer  Function 


4(S  +  2.5-) 

_ (S  +  10) _ 

3  ( S  +  3.33) 


(S  +  10) 


2(S  +  5) 


(S  +  10) 


10 


(S  +  10) 


(S  +  4) 


(S  +  2) 


(S  +  1) 


256 


(S2  +  22. 4S  +  256) 


144 


(S2  +  16. 8S  +  144) 


81 


(S2  +  12. 6S  +  81) 


36 


(S2  +  8.4S  +  36) 


16 


(S2  +  5.6S  +  16) 


512 


(S2  +  29.76S  +  256) (S2  +  12.16S  +  256) 
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Note:  0  -  indicates  control  system 


Appendix  B 


Stability  Derivatives  and  Aircraft  Transfer  Functions 

Data  from  the  Calspan,  NT-33  variable  stability 

aircraft  (Ref  2  ) ,  augmented  to  simulate  five  different 

aircraft  dynamics,  were  used  in  this  study.  This  appendix 

lists  the  equivalent  body  axis  stability  derivatives  of  the 

augmented  aircraft  (Table  Bl).  These  stability  derivatives 

were  estimated  during  a  prior  AFFDL  contract.  The  aircraft 

transfer  functions  —  ,  and  -y-  are  developed  using  these 

e  e  e 

derivatives.  The  subscript  used  with  the  transfer  functions 
denotes  the  aircraft  number.  Both  the  phugoid  and  the  short 
period  modes  are  included  in  the  first  set  of  equations.  The 
short  period  approximation  is  then  utilized  for  developing 
the  short  period  equations.  The  phugoid  and  short  period 
damping  ratios  and  natural  frequencies  are  evaluated  from  the 
resulting  equations  and  tabulated.  The  transfer  function 


and 


is  developed  as  follows: 

Flight  path  angle  y  is  given  by: 

y  =  0  -  a 

dividing  throughout  by  «  yields, 

JL  =  .9,  «  =  JL  -  JL 

6  «  6  6 
e  e  e  e 

Therefore  the  difference  of  the  transfer  functions 
is  taken  to  obtain  . 


(Bl) 


(B2) 

e 

6 
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Table :B1  NT-33  Aircraft  Dimensional  Stability  Derivatives 
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33685  0.33685  0.33685  0.33685  0.33685 


Substituting  the  values  of  stability  derivatives  and 
evaluating  the  determinants  for  each  aircraft  yields: 


e_  =  0. 336847S2  +  0.263892S  +  0.0195185 _ 

6e  1  S4  +  1.551S3  +  1.1353S2  +  0.072629S  +  0.018372 

e_  =  0 . 336847S2  +  0.262669S  +  0.0198526 _ 

6e  2  S4  +  2.671S3  +  5.4431S2  +  0.34444S  +  0.15437 

_e  =  0. 336847S2  -f  0.262025S  +  0.0198292 _ 

6e  3  S4  +  1.141S3  +  4.8991S2  +  0.29516S  +  0.18441 
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e_ 

s 

e 


4 


0.336847S2 
S4  +  4.211S 


+  0. 262519S  +  0.0198364 _ 

3  +  15.354S2  +  0.95023S  +  0.4887 


0. 336847S  +  0.258435S  +  0.0197128 

- >r - n - - — - 

S  +  4.211SJ  +  15.354S  +  0.95023S  +  0.4887 


and 


JL  =  S3  +  69.9339S2  +  4.11992S  +  2.666 _ 

6e  1  205(S4  +  1.551S3  +  1.1353S2  +  0.072429S  +  0.018372) 


a 

1.1SJ  + 

71.1108S^ 

4  4 .25147S 

4  2.77408 

6 

e 

2 

205 (S4  4 

2.671S3  4 

5.4431S  4 

0.3444S  4 

0.15437) 

a 

1.1S3  + 

69.4168S2 

4  4.1833S  4 

2.77408 

6 

e 

3 

205 (S4  4 

1.141S3  4 

4.8991S2  4 

0.29516S 

4  0.18441) 

a 

1.1S3  4 

69.4168S2 

4  4.32495S 

4  2.77408 

6 

e 

4  ~ 

205(S4  4 

4.291S3  4 

4.2149S2  4 

0.2 906 7S 

4  0.054619) 

a 

1.1S3  4 

72.6728S2 

4  4.31433S 

4  2.77408 

6 

e 

5 

205(S4  4 

4.211S3  4 

15.354S2  4 

0.95023S 

4  0.4887) 

The  transfer 


relationship 


fundtions,  , 
e 

£  =  _x_  are: 

~6 

e 


obtained  by  using  the 


JL  I  =  -S3  -  0.880265S2  +  49.97794S  n-  1.335243 _ 

|l  205(S4  +  1.551S3  +  1.1353S2  +  0.072429S  +  0.018372) 
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j__  =  -1.153  -  2.057165S2  +  49.59568S  +  1.2957 _ 

6e  2  205(S4  +  2.671S3  +  5.4431S2  +  0.3444S  +  0.15437) 


_y_  I  =  -1.1S3  -  0. 363165S2  +  49 . 53183S  +  1.291 

<5  I  o  /  a  ^  -i  r»  3  r  r\  o  O  c  -i  C  n 


e|  3  205(S  +  1.141S  +  4.899S  +  0.29516S  +  0.18441) 


_X_  =  -l.l?3  -  3.88317S2  +  49.491.445S  +  1.2924 _ 

6e  4  205(S4  +  4.291S3  +  4.2149S2  +  0.29067S  +  0.05462) 


JL  =  -1.1S3  -  3.5192S2  +  48 . 66485S  +  1.267 _ 

6e  5  205(S4  +  4.21153  +  15.354S2  +  0.95023S  +  0.4.887) 


The  perturbation  equations  using  the  short  period 
approximation  are: 


(S  -  Z  ) 
w 


-M  u 
wuo 


S(S  -  M  ) 


cx  0 

The  transfer  functions  —  and  —  are  given  by: 


Mt  S(S  -  M  ) 

e  ^ 
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-M  u 
w  o 


A  sp 


(B9) 


where , 


A  Sp 


(s  -  z)  -s 

w 

-Mu  S(S  -  M  ) 

wo  '  q 


(BIO) 


Substituting  the  values  of  stability  derivatives  the 
following  aircraft  short  period  transfer  functions  are 
obtained . 


0.33685S  +  0.25 _ 

S(S2  +  1.51S  +  1.04587) 


e_  =  0.33685S  +  0.25 _ 

6e  2  S(S2  +  2.63S  +  5.30673) 


0.33685S  +  0.2482 
S(S2  +  lvlS  +  4.8253) 


0.33685S  +  0.2487 _ 

S(S2  +  4.25S  +  4.012) 


e_  =  0.33685S  +  0.24463 
6e  5  S(S2  +  4.17S  +  15.15:5) 
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r 

. 

| 

a 

0.00537S  +  0.341 

6 

e 

1  (S2  +  1.513S+  1.04587) 

a 

0.00537S  +  0.3467 

6 

e 

2  (S2  +  2.63S  +  5.3067) 

a 

0.005  3  7S  +  0.3384 

6e 

3  (S2  +  1.1S  +  4.8253) 

a 

0.00537S  +  0.3556 

6 

e 

4  (S2  +  4.25S  +  4.012) 

a 

0.00537R  +  0.3543 

6 

e 

5  (S2  +  4.17S  +  15.155) 

The  transfer  functions  are  obtained  as  shown  earlier. 

e 

JL 

-0.00537S2  -  0.0041S  +  0.25 

6 

e 

1  S(S2  +  1.51S  +  1.04587) 

JL 

-0.00537S2  -  0.00982S  +  0.2489 

6 

e 

2  S(S2  +  2.63S  +  5.3067) 

JL 

-O.OQ537S2  -  O.OQ1556S  +  0.258 

6 

e 

3  S(S2  +  1.1S  +  4.8253) 

118 


-O.OQ5373S2  -  0.0187S  +  0.2487 
S(S2  +  4.25S  +  4.012) 


-O.OQ537S2  -  0.01744S  +  0.245 
S(S2  +  4.17S  +  15.155) 


Table  B2  Phugoid  and  Short  Period  and  K 


Aircraft 

No. 

Phugoid 

Short 

Period 

0) 

n 

5 

0) 

n 

K 

1 

0.132 

0.166 

1.026 

0.735 

2 

0.171 

0.148 

2.298 

0.57 

3 

0.196 

0.135 

2.19 

0.248 

4 

0.118 

0.25 

1.988 

1.064 

5 

0.18 

0.15 

3.886 

0.535 

Table  B3  Short  Period  u>n  And  £,' 


Aircraft 

No. 

Short  Period 
Approximation 

ID  „ 
n 

1 

1.023 

0.738 

2 

2.304 

0.571 

3 

2.197 

0.25 

4 

2.003 

1.061 

3.89 


0.536 


r 

I 

Appendix  C 

Method  of  Least  Squares  Curve  Fit 

This  appendix  presents  the  mathematical  development  of 
the  method  of  linear  least  squares,  used  for  the  analysis  of 
the  open-loop  frequency  response  data.  The  method  is  useful 
in  matching  an  arbitrary  set  of  data  points  with  a  known  set 
of  data  points.  The  mathematical  phylosophy  behind  the  method 
is  to  find  a  set  of  constants,  such  that  the  square  of  the 
error  between  a  known  set  of  data  and  a  linear  combination  of 
an  arbitrary  set  of  data  is  minimized.  The  known  set  of 
data  are  the  actual  pilot  ratings.  The  arbitrary  set  of  data 
points  are  the  phase  angles  and  slopes  of  the  phase  angles. 

The  method,  therefore,  is  applied  to  find  a  pilot  rating  ex¬ 
pression  as  a  linear  combination  of  the  phase  angle  and  the 
slope  of  the  phase  curve. 

Consider  a  set  of  three  data  points  u. ,  x.  and  y  whose 

11  i 

linear  combination  is  equal  to  PR.  Expressing  it  mathemati¬ 
cally  we  get: 

(PR)j  =  a,ui  +  alxi  +  a,y.  (Cl) 

where  a,,az  and  a3  are  constants.  For  our  problem  (PR)^  are 

the  estimated  pilot  ratings,  u^ ,  the  phase  angles,  and  x^,  y^ 

are  the  slopes.  If  7.^  represents  the  actual  pilot  ratings  , 
then  the  error  in  the  rating  expression  Ej^  is  given  by 
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Ei  =  Z.  -  (PR). 


Or 


and 


Ei  =  Zi  -  ( a , U±  +  a2Xi  +  a3yi) 


n  2  n 


ii/i  *  '  (°lUi  +  °2Xi  +  °3yl>}  - 


(C2) 


(C3) 


(C4) 


To  find  a,, a 2, a 5  such  that  eE^  is  a  minimum. 


3  i  =  l 


n 

5  Z[Z^  -  (c^u.  +  «2xi  +  “3y.)](-ui)  -  0  (C5) 


3  E 


2  n 


9  a 


=  £  2[Zi  -  (aiui  +  «2xi  +  “3yi)  K^)  =  O  (C6 ) 


2  i=l 


3E' 


n 


3a3  =  i!i  2[Zi  "  (“iUi  +  a^i  +  “3ui)](_yi)  =  0  (C7) 


Dividing  equations  C5,C6,C7  by  2  and  rearranging  yields, 

o 

XZiUi  =  (E^  )a1  +  (lu.x.)a2  +  (luiyi)a3  (C8) 

2 

lZixi  =  (Euixi)«1  +  (Ex.^  )a2  +  (Exiyi)a3  (C9) 


EZ 


iYi  =  <ruiyi>“i  +  (2xiyi)a2  +  >  3 


(CIO) 


Equations  C8,  C9  and  CIO  are  three  linear  algebraic  equa¬ 
tions  in  unknowns  a,,aj  and  a,.  Writing  the  equations  in  matrix 
format  yields: 
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Solving  for  a1}aX}a,  yields 


-1 


r  2 

ZU^  zu^x^ 

zu.  y . 

iJi 

—  *■ 

zZ.u. 

l  l 

- 

“  1 

„  r  2 

Zu.x.  Zx. 

11  1 

zx.  y . 

i7i 

zZ.x. 
i  l 

= 

zu^i  zxiy. 

^i_ 

lZl*l 

Equation  C12is  used  for  the 

evaluation 

Of  a, 

a  2  and  a  3 

(Cl  3 


develop  pilot  rating  expressions  described  in  chapter  III  of 


the  report.  It  can  be  modified  for  cases  where  there  are 
only  two  unknown  constants  a,  and  a2  by  deleting  the  third 
row  and  column. 


Appendix  D  Analog  Computer  Simulation  Diagrams 


This  appendix  contains  the  details  of  the  analog  computer 
circuitry  used  in  the  simulation  of  pilot-aircraft  dynamics 
and  the  control  system  dynamics.  Figure  D1  is  a  block  diagram 
description  of  the  system  simulated  as  a  closed-loop  pitch 
angle  tracking  task.  Figures  D2  and  D3  are  the  block  diagrams 
showing  closed-loop  flight  path  angle  tracking  task  without 
the  pitch  feedback  and  with  pitch  feedback  respectively.  For 
simplicity,  each  block  in  the  diagrams  are  discussed 
individually. 


Simulation  of  Aircraft  Short  Period  Dynamics 

The  aircraft  short  period  equations  of  motion  are  given 


by : 


a  -  7W«  -  6  = 


u 


-M  u  a  +  6  —  M  0  =  M.  S 
w  o  q  c>e  e 


Rearranging  these  equations  yields  the  following. 


« 

L  <*  +  0  +  '  6 

w  u  e 

o 


0=M0+Mua  +  M_  6 
q  w  o  6  e 


Y  =  0  -  a 


(Dl) 

(D2) 

(D3) 

(D4) 

(D5) 


Fig.  D1  Closed-Loop  Pitch  Tracking  Task 


Fig.  D3  Closed-Loop  y-Tracking  Task  With  '0'  Loop  Closed 


Eqns.  D3 ,  D4  and  D5  are  utilized  for  the  simulation.  By 
using  appropriate  values  of  the  dimensional  stability 


(Z  ,  M  .  M  )  and  control 
w  q  w 


u 


derivatives,  any  one  of 


o  e 

the  five  aircraft  configurations  can  be  easily  formed.  The 
analog  circuit  diagram  is  shown  in  Figure  D4. 


Simulation  of  Control  System 

Since  the  control  systems  varied  from  a  first  order  lead/ 
lag  to  a  first-order  and  second-order  lag  systems,  three  sep¬ 
arate  patching  areas  were  used,  one  for  each  type  of  control 
system.  The  three  categories  are  described  separately.  The 
fourth-order  lag  control  system  (c/S  13)  was  not  simulated. 

First-Order  Lag  Control  System 

A  first  order  lag  transfer  function  is  given  by 


v(s)  _  a 
x(s)  s  +  a 


(D6) 


where  x  and  y  are  the  input  and  output.  Rewriting  equation 
in  a  differential  equation  form  yields, 

y  =  -ay  +  ax  (D7) 

Equation  D7  can  be  simulated  as  shown  in  Figure  1)5.  The  value 
of  'a’  can  be  adjusted  to  simulate  any  desired  first-order  lag 
control  system. 
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Fig.  D4  Aircraft  Short  Period  Analog  Simulation  Diagram 


Fig.  D5  First-Order  Lag  Control  System  Simulation  Diagram 

126 


Second-Order  lag  Control  System 


A  second-order  lag  transfer  function  is  given  by, 

y(s)  _  _ b _ 

x(s)  S  +  as  +  b 

The  above  equation  can  be  written  as 


or 


ys/-  4  ays  +  by  =  bx 


_v  =  ,  by  _  bx 


(D8) 


(D9) 


(DIO) 


The  simulation  diagram  of  equation  DIOis  shown  in  Figure  D6 , 
where  ’a'  and  ’ b '  can  be  set  to  any  desired  value  to  simulate 
the  required  second-order  lag  system. 


Lead/Lag  Control  Svstem 

The  transfer  function  of  a  lead/lag  control  system  can 
be  written  as, 

v(s)  =  as  4  b  (Dll) 

x(s)  s  4  b 

Dividing  the  numerator  by  the  denominator  gives, 

xisl  ,  a  .  i>0  .1)  (D12) 

x(s)  (s  4  b) 

Rearranging  this  equation  yields 

y(s)  =  ax(s)  -  b(g  V)x(s)  (D3  3 ) 

(s  4  b) 
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Second  Order  Lag  Control  System  Simulation  Diagram 


Lead/Lag  Control  System  Simulation  Diagram 


The  simulation  diagram  of  eqn.  D13is  shown  in  Figure  D7. 
Appropriate  lead/lag  control  systems  can  be  simulated  by 
selecting  the  values  of  'a'  and  ’ b * . 


Simulation  of  Pilot  Model 

The  pilot  model  used  for  simulation  has  a  0.3  sec  delay, 
a  Dilot  gain  Kp,  and  a  fist-order  lead  circuit  with  pilot  lead 

time  constant  T^.  A  first  order  Pade  approximation  was  used 
for  the  time  delay  in  the  model.  The  pilot  transfer  function 


thus  reduces  to 


Y  (s)  =  Kp(TLs  +  1 )e 


-,s>  -yy  *  i><*  -  zm 

(s  +  2/t) 


For  ease  of  programming,  the  pilot  model  is  decomposed  into 
the  block  diagram  shown  in  Figure  D8  . 


Figure  D  8  Analog  Pilot  Model  Block  Diagram 


The  equation  for  the  first  block  is 


The  analog  computer  circuit  for  block  I  in  Figure  D  9  is 
shown  in  Figure  D9  . 


K 

P 

-  ei 

- 

o 

Figure  D9  Circuit  Diagram  For  Equation  D15 
The  transfer  function  of  the  second  block  (Figure  DIO)  is 


e^  s  -  2/t 
e2  s  +  2/t 


(D16) 


With  t  chosen  as  0.3  sec.,  Eqn«Dl6  maybe  rewritten  as 


6.667,  * 

-e-j  =  -e2  +  - (e2  +  e^) 


(D17 ) 


The  computer  simulation  diagram  of  Eqn.  D17  is  shown  in  Fig.  DIO. 


The  transfer  function  for  the  third  block  in  Figure  DIO  is 


^  =  tls  +  1  (D18) 

Eqn.  D18  shows  that  the  output  e^  consists  of  a  sum  of  the  in¬ 
put  and  its  derivative.  Simulation  of  derivatives  of  inputs 
is  not  desirable  in  analog  computer  programming.  Since  the 
value  of  used  was  about  0.5  sec.,  Eqn.  D18  was  modified  to 

include  a  pole  at  -20.  A  pole  at  -20  will  not  cause  any  sig¬ 
nificant  changes  and  the  errors  are  negligibly  small.  There¬ 
fore  Eqn.  D18  is  modified  as  shown 

e  a  (T.s  +  1)  ( 20T,  -  1) 

—  =  -  =  T. - (D19 ) 

e3  (s  +  20)  L  (s  +  20) 

The  simulation  diagram  of  Eqn.  D19  is  shown  in  Figure  Dll. 


Figure  Dll  Circuit  Diagram  For  Equation  D19 
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Proportional  Plus  Integral  Control 


The  proportional  plus  integral  control  was  used  in  the 
two-loop  closure  study.  In  Figure  D3,  with  the  inner-loop 
closed,  the  system  reduces  to  a  Type-o  System  and  therefore 
generates  a  steady-state  error  for  a  step  -y-command.  The 
proportional  plus  integral  control  makes  the  system  Type  1, 
thus  making  it  possible  for  the  output  to  follow  the  input  and 
driving  e^,  the  error  in  '  y',  to  zero. 

The  transfer  function  used  for  a  proportional  plus  integral 
control  is 

-  1  ♦  -f-  (D20 ) 

x(s)  s 

where  'a'  is  a  constant.  The  simulation  diagram  of  Eqn.D20  is 
shown  in  Figure  D12. 


Figure  D12  Circuit  Diagram  For  Equation  D20 
Simulation  of  ISE 


The  ISE  (integral  of  error  squared)  was  used  as  the  per¬ 
formance  criterion  for  adjusting  the  value  of  any  desired 
parameter.  The  simulation  diagram  for  ISE  is  shown  in  Figure 
D13.  Quarter  square  multipliers  were  used  for  squaring  the 


132 


Figure  D13  Quarter  Square  Multiplier  Simulation  Diagram 

The  analog  simulation  studies  were  thus  conducted  by 
patching  the  diagrams  discussed  above  on  one  patching  board. 
The  appropriate  blocks  were  cascaded  to  simulate  systems  shown 
in  Figures  D1 ,  D2  and  D3  . 
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